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A two-microelectrode voltage clamp and optical measurements of membrane potential changes at the transverse 
tubular system (TTS) were used to characterize delayed rectifier K currents (IK V ) in murine muscle fibers stained 
with the potentiometric dye di-8-ANEPPS. In intact fibers, IK V displays the canonical hallmarks of K v channels: 
voltage-dependent delayed activation and decay in time. The voltage dependence of the peak conductance (gK v ) 
was only accounted for by double Boltzmann fits, suggesting at least two channel contributions to IK V . Osmotically 
treated fibers showed significant disconnection of the TTS and displayed smaller IK V , but with similar voltage de- 
pendence and time decays to intact fibers. This suggests that inactivation may be responsible for most of the decay 
in IK V records. A two-channel model that faithfully simulates IK V records in osmotically treated fibers comprises a 
low threshold and steeply voltage-dependent channel (channel A), which contributes ~31% of gK v , and a more 
abundant high threshold channel (channel B), with shallower voltage dependence. Significant expression of the 
IKyl.4 and IK V 3.4 channels was demonstrated by immunoblotting. Rectangular depolarizing pulses elicited step- 
like di-8-ANEPPS transients in intact fibers rendered electrically passive. In contrast, activation of IK V resulted in 
time- and voltage-dependent attenuations in optical transients that coincided in time with the peaks of IK V records. 
Normalized peak attenuations showed the same voltage dependence as peak IK V plots. A radial cable model includ- 
ing channels A and B and K diffusion in the TTS was used to simulate IK V and average TTS voltage changes. Model 
predictions and experimental data were compared to determine what fraction of gK v in the TTS accounted si- 
multaneously for the electrical and optical data. Best predictions suggest that K v channels are approximately 
equally distributed in the sarcolemma and TTS membranes; under these conditions, >70% of IK V arises from 
the TTS. 



INTRODUCTION 

Voltage-dependent delayed rectifier K channels (K v ) are 
known to play a crucial role in skeletal muscle physio- 
logy; they are responsible for the downstroke phase of 
the action potential (AP) that rapidly reestablishes the 
resting membrane potential after the opening of Na 
channels. The overall properties of K v currents have 
been mostly studied in muscle fibers from the frog 
(Adrian et al., 1970; Adrian and Marshall, 1976) and the 
rat (Duval and Leoty, 1980; Pappone, 1980; Beam and 
Donaldson, 1983a,b) and to a much lesser extent in fibers 
from the mouse (Brinkmeier et al., 1991; Hocherman 
and Bezanilla, 1996). The studies in mouse fibers have 
limitations derived from the fact that they have been per- 
formed using several configurations of the patch-clamp 
technique. For example, when on-cell or excised patch 
configurations were used (Hocherman and Bezanilla, 
1996), no information was obtained about K channels 
potentially located in the transverse tubular system (TTS) 
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membranes or about the ensemble properties of currents 
from the entire muscle cell. Alternatively, attempts to eval- 
uate the properties of K v currents (IK V ) using the whole- 
cell patch-clamp configuration (Brinkmeier et al., 1991) 
suffer from technical limitations possibly related to the 
large magnitude of the currents. Consequently, a more 
detailed characterization of IK V in the mouse is timely. 
The application of the two-microelectrode voltage-clamp 
technique in short fibers from the foot muscles of the 
mouse (flexor digitorum brevis [FDB] or interosseous 
muscles) is currently accepted as the most adequate 
approach to investigate the electrophysiological prop- 
erties of muscle fibers without the aforementioned 
limitations (Friedrich et al., 1999; Ursu et al., 2004; 
DiFranco etal., 2011a; Fu etal., 2011). 

It is generally postulated that IK V in adult mammalian 
muscle fibers display decaying phases that result from 
channel inactivation and/ or K accumulation in the 
lumen of the TTS, indirecdy implying that a fraction of 
K v channels may be located in the TTS. Thus, though 
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the presence of IK V contributions arising from both the 
TTS and surface membranes has been suggested for 
rat skeletal muscle (Duval and Leoty, 1980; Beam and 
Donaldson, 1983a), no specific information regarding the 
K v channel distribution is available in the literature. 

The identification of K v channels in skeletal muscle 
has been undertaken mostly using molecular biology 
and biochemical approaches. Using Northern blotting 
analysis, several types of K v channels have been identi- 
fied in adult mice, including members of the Shaker 
(e.g., K v l.l, K v 1.4, Kyl.5, and K v 1.7) and Shaw (K v 3.1 
and K v 3.4) subfamilies (Lesage et al., 1992; Kalman 
et al., 1998; Vullhorst et al., 1998) and members of the 
slowly activating and inactivating K v subfamily (K v 7.2, 
K v 7.3, and K v 7.4; Iannotti et al., 2010). Nevertheless, 
only K v 3.4 and K v 1.5 have been reported to be ex- 
pressed (as proteins) in rat and human muscles (Abbott 
et al., 2001; Bielanska et al., 2009). Interestingly, re- 
cent reviews about ionic channel genes expressed in 
skeletal muscle membranes suggest that only K v 1.4, 
K v 3.4, and K v 7.4 may be functionally important in this 
tissue, but no evidence supporting this statement is 
given (Jurkat-Rott et al., 2006; Kristensen and Juel, 
2010). Although the currents carried by K v isoforms 
expressed in heterologous systems have been studied 
(Po et al., 1993; Abbott et al., 2001), limitations of this 
approach weaken the implications for native K v cur- 
rents in adult muscle fibers. For example, it is well 
known that K v channels are assembled in vivo from 
more than one a subunit isoform (Ruppersberg et al., 
1990; Po et al., 1993) and that tetramers are regulated 
by accessory subunits (Abbott et al., 2001; Pongs and 
Schwarz, 2010). To our knowledge, there are no pub- 
lished attempts comparing properties of IK V recorded 
from adult muscle fibers and those from heterologous 
expression systems, but recent studies about the pro- 
perties of Nayl-4 channels suggest that results in vitro are 
not readily applicable to in vivo conditions (DiFranco 
and Vergara, 2011; Fu et al., 2011). Considering the 
crucial role that IK V play in skeletal muscle fibers, the 
scarcity of information is surprising, and the present 
work aims to fill the gap. 

We first determine the kinetic and voltage depen- 
dence of IK V and decide whether the currents are trans- 
ported through one or more K v channels. For this, we 
compare the properties of IK V in intact and osmotically 
treated (detubulated) muscle fibers. We also provide evi- 
dence about the molecular identity of K v channels that 
are expressed in the muscles from which the fibers used 
for electrophysiological measurements are isolated. Fi- 
nally, we address the question of the relative distribu- 
tion of K v channels between the TTS and surface 
membranes. To this end, we simultaneously record IK V 
and fluorescence transients from fibers stained with 
the potentiometric indicator di-8-ANNEPS and, as re- 
cendy reported for other conductances (DiFranco and 



Vergara, 2011; DiFranco et al., 2011a), use radial cable 
model simulations of the membrane potential changes 
in the TTS to predict the relative distribution of K v chan- 
nels between the surface and TTS membranes within 
a narrow range. 

MATERIALS AND METHODS 

Biological preparation 

Animal handling followed die guidelines laid down by the University 
of California, Los Angeles, Animal Care Committee. FDB and in- 
terossei muscles from C57BL 3-4-mo-old mice were used. Isolated 
fibers were obtained by enzymatic dissociation as previously de- 
scribed (Woods et al., 2004; DiFranco et al., 2011a). Two groups 
of fibers were used in electrophysiological experiments: control 
fibers, which were maintained intact after isolation (n = 44), and 
osmotically shocked fibers, which were selected after formamide- 
based osmotic treatment aimed to disconnect the TTS from the 
surface membrane (n = 40). From this latter group, only those 
displaying a specific capacitance of <2.5 uF/ cm 2 were further ana- 
lyzed (n = 12). The geometrical and electrical parameters of both 
intact and osmotically treated fibers are shown in Table 1 . 

Solutions 

The composition of the solutions (in mM) was as follows. Tyrode: 
156 NaCl, 10 MOPS, 2 CaCl 2 , 10 dextrose, 1 MgCl 2 , and 4 KC1, 
pH adjusted with NaOH. NMG-Tyrode: 150 NMG, 5 RbCl, 10 MOPS, 
2 CaCl 2 , 10 dextrose, 1 MgCl 2) 4 KC1, and 2 x 10~ 4 tetrodotoxin 
(TTX), pH adjusted with HC1. TEA-Tyrode: 145 TEA-OH, 10 MOPS, 
10 CsOH, 2 Ca(OH) 2 , 1 Mg(OH) 2 , 5 dextrose, 2 x 10~ 2 nifedi- 
pine, 0.4 9-anthacene carboxylic acid (9-ACA) , and 2 x 10~ 4 TTX, 
pH adjusted with HC1. Formamide-Tyrode: Tyrode added with 2 M 
formamide. K internal solution: 160 KOH, 20 MOPS, 50 EGTA, 
5 ATP-Mg, 5 Na 2 -creatine phosphate, and 5 of reduced glutathi- 
one, pH adjusted with aspartic acid. The high EGTA concentra- 
tion in the internal solution was used to arrest fiber contraction 
to avoid movement of artifacts in the optical records. The osmo- 
larity of the aforementioned solutions was 300 ± 5 mosmol/kg 
H 2 0. Tris-buffered saline (TBS): 10 Tris-Cl and 150 NaCl. Lysis 
buffer: 3% SDS, 115 sucrose, and 66 Tris-Cl. Fractionation buf- 
fer: 150 KC1, 5 MgS0 4 , 20 MOPS, and 10 EGTA-K. Microsome 
storage buffer: 20 MOPS-Tris and 300 sucrose. Stock solutions 
for 9-ACA (500 mM), nifedipine (50 mM), and isradipine (10 mM) 
were prepared in DMSO. All solutions were adjusted to pH 7.4. 
All chemicals were purchased from Sigma-Aldrich. 

Osmotic shock treatment 

The method was based on the use of formamide (del Castillo and 
Escalona de Motta, 1978) and modified from that described else- 
where for enzymatically dissociated mouse FDB fibers (Lueck et al., 
2010; DiFranco and Vergara, 2011). In brief, dissociated fibers 
were transferred to a 50-ml glass beaker containing 0.5 ml Tyrode, 



TABLE 1 

Geometrical and electrical properties of muscle fibers 



Property 


Intact fibers 


Osmotically treated 




(« = 44) 


detubulated fibers ( n= 11) 


Radius (p) 


254 ± 3.7 


24.8 + 3.0 


Length (um) 


492 ± 80 


501 ± 57 


Capacitance (LiF/cm^) 


5.00 ± 0.49 


2.3 + 0.31 a 



Values correspond to the mean ± SD. 
Significance at P < 0.05. 
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and 10 ml formamide-Tyrode was added. After 10 min, the solu- 
tion volume was reduced again to M).5 ml, and 40 ml Tyrode was 
quickly added. Approximately 10 min after the osmotic shock, the 
fibers were stained with di-8-ANEPPS as described previously 
(DiFranco and Vergara, 2011). The criteria used to evaluate the 
extent of detubulation are described in Results (see Osmotic 
shock treatment affords a degree . . .). 

Electrophysiology 

Experiments were performed under voltage-clamp conditions 
using a two-microelectrode high voltage amplifier (TEV-200A; 
Dagan Corporation) as previously described (Woods et al., 2004, 
2005; DiFranco et al., 2005, 2011a). To improve the frequency re- 
sponse of the voltage clamp, electrodes were drawn to the largest 
tip size compatible with the viability of the fibers and had resis- 
tances in the range of 6 to 10 Mfi when filled with internal solu- 
tion. Experiments with detubulated fibers required the use of 
electrodes in the lowest range of resistance to ensure voltage-clamp 
stability. Command voltage pulses were digitally filtered at 50 kHz. 
Fibers were impaled under current-clamp conditions and, after an 
^20-min equilibration period between the pipette solution and 
the myoplasm, were voltage clamped at a holding potential (V H ) of 
—90 mV. Fibers requiring >12 nA to maintain Vh were discarded. 
Two protocols were used to eliminate capacitive and linear ionic 
currents. (1) Current records elicited by 40-50-mV pulses was 
linearly scaled and subtracted from the rest of the records. (2) 
Current records elicited by any given pulse in TEA-Tyrode were 
subtracted one-to-one from those in which nonlinear currents 
were recorded. In contrast to the first protocol, the second proto- 
col subtracts not only linear components, but also nonlinear ca- 
pacitive currents. Linear membrane capacitance was measured 
in every fiber as previously described (DiFranco et al., 2011a) for 
both intact and osmotically treated fibers. Potassium currents were 
typically expressed in microamperes/ centimeters squared. In 
some cases, results were also expressly given in amperes/farads for 
comparison, but a general conversion to this unit can be readily 
made using the mean specific capacitance in Table 1 . 

The voltage dependence of the peak gK v was adjusted with either 
a single Boltzmann equation or by the following double function: 



For control fibers, the illumination disc was adjusted to ~25 pm 
in diameter and focused (using a microscope [UPLSAPO 
100XO; Olympus] and a 1.4 numerical aperture objective) at the 
center of the x, y, and z axes of the fiber. For detubulated fibers, 
the illumination spot was widened to a diameter ~25% larger 
than the fiber's width to include the periphery of the fibers. In either 
case, the illumination spot was centered at the site where the 
voltage microelectrode was impaled. Optical signals were low-pass 
filtered (2 kHz) single sweeps. Optical data are presented as nor- 
malized fluorescence changes (AF/F). 

Radial cable model simulations 

The properties of K currents and the impact that they have on mem- 
brane potential changes in the TTS were simulated using multistate 
models and the radial cable model, respectively. Details of model 
equations and integration methods are described in the Appendix. 

Total protein lysates and microsomal fractions 

To prepare total protein lysates from foot muscles (muscle lysate 
[ML] ) , FDB and interossei muscles from adult mice were minced 
with a razor blade and homogenized on ice using a tissue grinder 
(Duall; Kontes) in lysis buffer with one tablet/7 ml of a cocktail of 
protease inhibitors (Roche). The muscle homogenate was rotated 
at 4°C for 1 h and then centrifuged at 15,000 gior 15 min. The su- 
pernatant was collected and stored at — 80°C. For muscle micro- 
somal (MM) preparations, hind limb muscles from two mice were 
dissected and trimmed of connective tissue and fat. Muscles were 
finely minced using razor blades, mixed with ice-cold homogeni- 
zation buffer supplemented with protease inhibitor cocktail (4 ml/ g 
of wet tissue), and homogenized for 30 s (Bio-homogenizer 
M133/1281; Biospec Products Inc.). The homogenate was centri- 
fuged at 1,500 g-for .30 min at 4°C. The supernatant was saved, and 
the pellet was resuspended, homogenized, and centrifuged again 
at 1,500 gfor 30 min. The two supernatants were combined and 
centrifuged at 15,000 gfor 30 min at 4°C. The new supernatant was 
centrifuged at 125,000 gfor 60 min at 4°C. The final pellet was re- 
suspended in a small volume of storage buffer and stored at — 80°C. 
For both preparations, protein concentration was determined using 
the Quick Start Bradford Protein Assay (Bio-Rad Laboratories) . 



peak gK v (V) = A 



1- 



1 + e 



■(V-Vl)/kl 



) (1 



-(V-V2)/k2 



where A is the peak [gK v ] m ax and e is the fraction of conductance 
contributed by one channel. For statistical analysis, individual da- 
tasets obtained at various conditions were fitted, and the means 
for each parameter were comparatively evaluated. All experi- 
ments were performed at room temperature (20-22°C). 

Recording of di-8-ANEPPS transients 

The optical methodology has been described elsewhere (DiFranco 
and Vergara, 2011; DiFranco et al., 2005, 2011a). In brief, disso- 
ciated fibers stained with di-8-ANEPPS (Biotium; Hayward) were 
placed on coverslip-bottomed small Petri dishes sitting on the 
stage of an inverted microscope (IX-71; Olympus) equipped with 
a standard epifluorescence attachment, a cooled charge-coupled 
device camera (ST-402ME; Santa Barbara Instrument Group), 
and a photodetector consisting of a photodiode (UV-001; OSI 
Optoelectronics) connected to a patch-clamp amplifier (Axopatch 
2A; Molecular Devices) . For control fibers, only those displaying a 
sharp sarcomere banding and the distinctive pattern of TTS stain- 
ing by di-8-ANEPPS (DiFranco et al., 2005, 2007) were used for 
the experiments. For detubulated fibers, only those displaying 
predominant peripheral di-8-ANEPPS staining were selected. 



Western blotting 

Western blot analysis was performed in three total protein MLs, two 
MM preparations, and two hippocampus lysates (HLs; Delgado and 
O'Dell, 2005). Gels were loaded with 40, 32, and 35 pg of total pro- 
tein for ML, MM, and HL samples, respectively. Proteins were sepa- 
rated on 7% SDS-PAGE gels (100 V for 1 h) and transferred to 
polyvinylidene fluoride membranes (25 V overnight). Membranes 
were blocked for 1 h in blocking buffer (4% skimmed milk powder 
and 2% bovine serum albumin in TBS with 0.05% Tween [TBS-T] ) . 
Immunoblotting was performed overnight using the following anti- 
bodies diluted in TBS-T containing 4% nonfat milk: mouse mono- 
clonal anti-K v 1.4 (75-010, clone K13/31, 1:500 dilution), mouse 
monoclonal anti-K v 3.4 (75-112, 1:500 dilution), and rabbit poly- 
clonal anti-K v l .4 (APC-007, 1 :200 dilution) . The first two antibodies 
were obtained from the University of California, Davis/National In- 
stitutes of Health (NIH) Neuromab Facility, and the third one was 
obtained from Alomone Laboratories. After three washes for 3 min 
in TBS-T, membranes were exposed to horseradish peroxidase-con- 
jugated secondary antibody (goat anti-mouse or goat anti-rabbit 
antibody, as required; MP Biomedicals; 1:2,000 dilution in TBS-T) 
for 2 h. Enhanced chemiluminescence with Immun-Star horse- 
radish peroxidase substrate (Bio-Rad Laboratories) was used to 
develop all immunoblots. Images were collected using a charge- 
coupled device camera attached to a ChemiDoc chemilumines- 
cent detection system and using Quantity One software (Bio- 
Rad Laboratories). 
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Data acquisition and statistical analysis 

Voltage, current, and fluorescence records were filtered at 10, 
5, and 2 kHz, respectively, using multiple pole analogue Bessel 
filters. Data points were sampled every 30 ps using a data acqui- 
sition interface (PCI-6221; National Instruments) and custom 
software written in LabView (National Instruments). Unless other- 
wise stated, pooled data are expressed as means ± SEM. Signifi- 
cance was set at P < 0.05. The goodness of the fits comparing 
single and double Boltzmann functions was evaluated using the 
Akaike's Information Criterion (AIC) test (Origin Pro8; Origin- 
Lab Corporation) . 

RESULTS 

Delayed rectifier potassium currents (IK V ) 
in intact FDB muscle fibers 

Fig. 1 A shows the ionic currents recorded in an FDB 
fiber under conditions that eliminate the contribu- 
tions (by specific ion replacements and/ or the addi- 
tion of blockers in the external solution) from Na 
channels (Nayl.4; by replacing Na by NMG and add- 
ing 4 x 10 7 M TTX) and CI channels (C1C-1, by add- 
ing 400 pM 9-ACA) and reduce those from inward 
rectifier K channels (Kj R ; by adding 5 mM Rb) . The 
fiber was voltage clamped at a holding potential (Vh) 
of -90 mV (close to the calculated E K of -91 mV), 
and 50-ms rectangular voltage pulses (from 0 to 
200 mV every 10 mV) were applied. The current re- 
cords were corrected by subtracting scaled linear leak 
and capacitive currents obtained from a record in re- 
sponse to a 50-mV pulse. Negligible ionic currents 
were observed in response to hyperpolarizing pulses 
(not depicted) or to pulses smaller than 70 mV (Fig. 1 A, 
gray trace). For larger depolarizations, outward cur- 
rents with typical features of IK V were observed: cur- 
rents arise after a delay from the pulse onset, and the 
larger the pulse amplitude, the shorter the delay be- 
comes. A very early nonlinear capacitive current com- 
ponent (gating currents) precedes the delayed ionic 
component. At every depolarization, IK V grows to an 
apparent maximum with a rate of rise that increases 
with the pulse magnitude. For the largest pulse tested 
(200 mV), a peak IK V of 1,100 pA/ cm was reached in 
^3.7 ms (Fig. 1 A, orange trace) . For physiological de- 
polarizations (to membrane potentials of ^40 mV), 
similar to those reached at the peak of an AP, a peak 
IKy of ^460 pA/ cm 2 is attained in ^5 ms (Fig. 1 A, 
blue trace) . Within the timescale used in Fig. 1 A, only 
IK V records elicited by relatively large depolarizations 
(>100 mV) already show evidence of decay during 
the pulse. 

The peak values of the IKy records shown in Fig. 1 A 
(peak IKy) are plotted as a function of membrane po- 
tential in Fig. 1 B (closed circles). It can be observed 
that peak IKy increases almost linearly with membrane 
depolarizations for the entire range of voltages ex- 
plored. Assuming an ohmic behavior for open channels 



(Hodgkin and Huxley, 1952a) , we calculated the peak K 
conductance (peak gK v ) as 

peak g K v =«l. 

The resulting peak gK v values are plotted as a function 
of the membrane potential in Fig. 1 B (open circles). 
Notwithstanding the ample voltage range explored, 
peak gK v did not reach a plateau; a maximum value of 
^5 mS/ cm was obtained in this fiber for a 200-mV 
pulse. It is also apparent that the voltage dependence 
of the data points do not seem to conform to a single 
sigmoid distribution, but instead there seems to be a 
noticeable inflection point at ^25 mV, potentially sug- 
gesting the participation of more than one channel's 
population. The voltage dependence of the mean 
peak IKy and gK v obtained from 16 experiments in 14 
fibers, is shown in Fig. 1 (C and D, respectively). The 
mean membrane potential at which sizable currents 
were detected is —20 mV, and the mean maximum 
peak IK V at 1 10 mV is 990 ± 193 pA/cm 2 (mean±SD), 
which corresponds to ^191 ± 37 A/F. The population 
data in Fig. 1 C display similar features as those from 
the individual fiber in Fig. 1 B. The peak conductance 
plot, calculated from mean data in Fig. 1 C, is shown in 
Fig. 1 D (circles). As observed in Fig. 1 B, a plateau is 
not reached within the voltage range explored. The 
maximal peak gK v ((gK v ) max ) obtained from the mean 
plot is 5 ± 1 mS/ cm 2 (mean ± SD) , which corresponds 
to ^1 ± 0.2 mS/pF. As suggested from data in Fig. 1 B, 
the voltage dependence of the mean data (open cir- 
cles) cannot be correctly fitted to a single Boltzmann 
dependence (Fig. 1 A, green trace). In contrast, a dou- 
ble Boltzmann curve (Fig. 1 D, black trace) adequately 
predicts the data; the AIC test indicates that this model 
function is >2 x 10 8 times better in predicting this 
mean dataset than the single function. This feature is 
suggestive that there are at least two populations of K 
channels contributing to the total IKy records in this 
preparation. One contribution seemingly arises from a 
minority population of low threshold channels with a 
half-activation voltage (V1/2) at approximately —10 mV 
(Boltzmann fit represented by the red trace in Fig. 1 D) , 
and the other arises from a larger population of higher 
threshold channels with V1/2 at 40 mV (Boltzmann fit 
represented by the blue trace in Fig. 1 D) . Interest- 
ingly, the two single Boltzmann curves in Fig. 1 D are 
sufficiently distinct to allow us to infer that the low thresh- 
old channel (red trace) is fully activated at approxima- 
tely — 5 mV, shows steep voltage dependence, and 
contributes 40% of the total conductance. In contrast, 
the high threshold channel (Fig. 1 A, blue trace) dis- 
plays shallower voltage dependence, is fully activated 
at a potential >80 mV, and contributes ^60% of the 
total current. 
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Figure 1 . Delayed rectifier K currents in FDB fibers. (A) IK V recorded in response to depolarizations from Vh to membrane potential 
values ranging from —30 to 110 mV in 20-mV steps (gray to orange traces). The NMG-Tyrode contained no Ca 2+ channel blockers. Lin- 
ear leak and capacitances were removed from each raw current record by subtracting a scaled version of the current record obtained 
in response to a 50-mV pulse. (B) Voltage dependence of the peak IK V (closed circles) and peak gK v (open circles) . For both plots, the 
circles are connected with straight lines. (C) Voltage dependence of the mean peak IK V obtained from 16 voltage families in 14 different 
fibers. (D) Voltage dependence of the mean peak gK v (circles), calculated from the data in C. The green trace is a single Boltzmann 
fit to the data with the following parameters (mean ± SD): A = 4.77 ± 0.9 mS/cm 2 ; VI = 18.2 ± 6.3 mV; Kl = 20.4 ± 1.5 mV. The black 
trace is a double Boltzmann fit to the data with the following parameters: A = 4.93 ± 0.9 mS/cm 2 ; e = 0.36 ± 0.06; VI = —7 ± 4.7 mV; Kl = 6.7 ± 
1.1 mV; V2 = 39 ± 6.6 mV; K2 = 19.9 ± 2.3 mV. The red and blue traces are the independent components of the double Boltzmann func- 
tion. Error bars represent SEM. 



IK V records from mouse FDB fibers show a marked decay 
during long depolarizations 

Although the timescale used in Fig. 1 was adequate to 
illustrate the activation characteristics of IK V , it was 
insufficient to study the voltage and time dependence 
of the decay processes. When a longer time window is 
used, as illustrated in Fig. 2 A, pronounced decays in 
the currents are observed. Furthermore, it is apparent 
from IKy records in Fig. 2 A that the decay process is 
nonmonotonic. Another striking feature of the family 
of currents in Fig. 2 A is the presence of prominent 
inward current tails at the end of the pulses, which 
rapidly decay toward the baseline. Because the fiber 
was repolarized to —90 mV, a value close to the calcu- 
lated E K , negligible current contributions through 
IK V channels would be expected. A possible source for 
these tail currents is, as suggested previously (Beam 
and Donaldson, 1983a), that K accumulation had oc- 
curred in the TTS and that some of the current con- 
tributions arise from this membrane compartment 
through K v channels (while they close) or through K IR 
channels. Because 5 mM Rb reduces the latter (Beam 



and Donaldson, 1983a) while not affecting IK V (un- 
published data), at least part of the tail currents may 
represent the closing of K v channels in the TTS. In 
addition, Ca 2+ current tails, through the prominent 
Cayl.l channels in the TTS, could represent an im- 
portant contribution to the observed currents at the 
end of the long pulses shown in Fig. 2 A, particularly 
after large depolarizations. It should be noted that this 
particular experiment was conducted with 2 mM Ca 2+ 
in the external solution, and no Ca 2+ channel blocker 
was used. In fact, slight distortions suggestive of Ca 2+ 
current contaminations can be observed as minor down- 
ward deflections in the IK V traces at large depolarizations 
(Fig. 2 A, deep blue and orange traces) . To remove the 
possible contributions of putative Ca 2+ currents and 
obtain IKy in complete isolation from other currents, 
we used known Ca 2+ blockers such as isradipine and 
nifedipine (Lamb and Walsh, 1987; Berjukowetal., 2000; 
DiFranco et al., 2011b). Fig. 2 B shows IKy records 
from the same fiber but in the presence of 0.5 jiM 
isradipine. It can be observed that a large fraction of the 
tail currents was blocked by isradipine. However, in 
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Figure 2. IK V decays and blocking effect of isradipine. (A) IK V records were obtained in response to 500-ms voltage pulses of 60-200-mV 
amplitude, with 20-mV increments. (B) IK V from the same fiber as in A, but after exposure to 0.5 uM isradipine. (C) Voltage dependence 
of peak IK V (closed circles) and peak gK v (open circles) determined from the data in A (black circles, control conditions) and B (red 
circles, isradipine) . The black solid lines are double Boltzmann fits to the control data with the following parameters: A = 4.44 mS/ cm 2 ; 
e = 0.33; VI = -3.2 mV; Kl = 6.2 mV; V2 = 44.7 mV; K2 = 17.4 mV. The red solid lines are double Boltzmann fits to the isradipine data 
with the following parameters: A = 2.68 mS/cm 2 ; e = 0.31; VI = -1.8 mV; Kl = 7.0 mV; V2 = 36.7 mV; K2 = 14.8 mV. (D) Mean peak gK v 
(pooled data) obtained under control conditions (black circles; same data as in Fig. 1 D) and in the presence of isradipine (red circles; 
n = 6) . The red line is a double Boltzmann fit to the isradipine data with the following parameters (mean ± SD) : A = 2.90 ± 0.7 mS/cm 2 (*) ; 
e = 0.28 ± 0.09 (*); VI = -1.7 ± 2.8 mV (*); Kl =7.6 ± 1.5 mV; V2 =40 ± 6 mV; K2 = 20.7 ± 3.5 mV. The * indicates statistical significance 
with respect to the control (P < 0.05). Error bars represent SEM. 



agreement with results obtained with other Ca 2+ channel 
blockers in other preparations (Grissmer et al., 1994; 
Zhang et al., 1997), Fig. 2 B shows that isradipine also 
has an important apparent blocking effect on IK V . In 
our case, ^40% of the maximal peak currentwas blocked 
by the drug. Importandy, the effects of isradipine on IK V 
records seem to be independent of the entry of Ca 2+ 
through Cayl-1 because IK V reaches a peak at times 
when activation of the Ca 2+ channels is negligible (un- 
published data). By comparing the records in Fig. 2 
(A and B), it is also noticeable that IK V decays to a 
larger extent in the presence of isradipine. This is an 
intriguing observation because the elimination of an 
inward current component (ICa) would have been ex- 
pected to have the opposite effect; we suggest that a 
reasonable explanation for this phenomenon is that 
the Ca 2+ entry through Cayl.l channels facilitates the 
activation of a significant outward K current through 
BK channels in the TTS (Latorre et al., 1983; Tricarico 
et al., 1997). The voltage dependence of the peak IK V 



and peak gK v before (black circles) and after exposure 
to 0.5 pM isradipine (red circles) is shown in Fig. 2 C. 
It can be seen that, for all voltages tested, isradipine 
reduces the peak gK v . This suggests this drug blocks 
the two populations of channels proposed above; in 
support of this, we found that, as in the case of control 
peak gK v , double Boltzmann equations were required 
to accurately fit the data obtained in the presence of 
isradipine (the AIC test indicates that this function is >9 x 
10 5 better than a single Boltzmann function; Fig. 2 C, 
red traces) . Nevertheless, the results obtained from a 
population of fibers (Fig. 2 D), which include data 
from 16 fibers in control conditions and 6 fibers in 
the presence of 0.5 pM isradipine, suggest that the pu- 
tative channels may be differentially affected by this 
drug. Thus, statistical analysis of the mean parameters 
of the Boltzmann fits to the peak gK v in Fig. 2 D shows 
that not only (gK v ) max , but also the percent contribu- 
tion of the low threshold channel are significantly re- 
duced by isradipine. These results imply that although 
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both channels are blocked by isradipine, the low thresh- 
old channel is more sensitive to the drug that the 
higher threshold one. 

Osmotic shock treatment affords a degree of TTS 
disconnection that allows for better measurements of IK V 
The possibility that two populations of K channels are 
contributing to IK V and the difficulty in eliminating ICa 
contaminations without affecting IK V complicate the 
analysis of this current in skeletal muscle fibers. For 
example, as shown for the Na channel (DiFranco and 
Vergara, 201 1 ) , current contributions from K v channels 
in the TTS are expected to have a different overall im- 
pact on IK V than those at the surface membrane. In ad- 
dition, though the Cayl.1 conductance is likely to be 
smaller than that of K v channels, its preferential loca- 
tion at the TTS membranes could be responsible for the 
distortion seen in IK V records. For these reasons and in 
the hope that a certain degree of disconnection of the 
TTS from the sarcolemma affords better isolation of 
IK V , as shown previously for the Na channel (DiFranco 
and Vergara, 2011; Fu et al., 2011), we submitted FDB 
fibers to a formamide osmotic shock. 

The extent of TTS disconnection was evaluated by 
two-photon laser-scanning microscopy (TPLSM) and 
electrophysiological methods. Fig. 3 A is a TPLSM image 
that illustrates the typical staining pattern of the sur- 
face and TTS membranes in a control FDB muscle 
stained with di-8-ANEPPS. As described in previous 
publications from our laboratory (DiFranco et al., 2005, 
2007, 2009), the TTS is imaged as paired bands of high 
fluorescence intensity (two T tubules per sarcomere) 
running perpendicular to the sarcolemma. These paired 
bands repeat regularly, with similar intensity, along the 
longitudinal axis of the fibers (consecutive pairs every 
2.9 um) . The inset in Fig. 3 A shows an expanded view 
of the surface and TTS staining that allows for the 
identification of the different distances separating the 
members of each pair of adjacent T tubules with respect 
to those separating consecutive T tubule pairs. Equiva- 
lent images from an FDB fiber stained with di-8-ANEPPS 
after being osmotically shocked are shown in Fig. 3 B. 
It is clear that, in contrast to Fig. 3 A, the TTS staining 
here is reduced to the periphery of the fiber, as indi- 
cated by the lack of "internal" staining. The staining 
pattern displayed in Fig. 3 B is consistent with the 
notion that osmotic treatment restricted the access of 
di-8-ANEPPS only to the periphery of the TTS; in other 
words, the bulk of the TTS is effectively disconnected 
from the surface membrane. However, TTS "stumps" 
of variable length can be readily observed up to ^2-3 um 
from the surface membrane. This is better appreciated 
in tangential sections as shown in the inset of Fig. 3 B, 
which shows that the pattern of TTS staining is pre- 
served at the periphery of the fiber. The pattern of 
di-8-ANNEPS staining of the TTS in TPLSM images of 
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Figure 3. TPLSM evaluation of TTS disruption by osmotic shock 
treatment. (A) TPLSM fluorescence image of an intact FDB mus- 
cle stained with di-8-ANEPPS. The inset is a threefold enlarge- 
ment of the area defined by the white rectangle. (B) TPLSM 
fluorescence image of an isolated fiber stained with di-8-ANEPPS 
after being submitted to a formamide-based osmotic shock. In the 
area delimited by the rectangle, the plane of the TPLSM section 
was tangential to the fiber periphery. This area is enlarged three- 
fold in the inset. 



control muscle fibers (as in Fig. 3 A) has been consis- 
tentiy established. In contrast, di-8-ANNEPS staining pat- 
terns like that illustrated in Fig. 3 B are observed in <15% 
of osmotically treated fibers; the rest of the treated fibers 
showed images of TTS staining with intermediate fea- 
tures between those of Fig. 3 (A and B) . Membrane ca- 
pacitance measurements provide support to the TPLSM 
data. The capacitance of a subpopulation of osmotically 
treated fibers selected for electrophysiological experi- 
ments was as low as 2.3 ± 0.31 pF/cm 2 (n= ll);however, 
the majority of treated fibers that were electrically tested 
had larger capacitances (unpublished data) . It is impor- 
tant to note that, other than differences in capacitance, 
untreated and osmotically treated fibers had indistin- 
guishable geometrical features (Table 1). 

Fig. 4 A and its inset show, at two timescales, IK V re- 
cords obtained from an osmotically shocked fiber in 
response to depolarizations (in 10-V steps) from —90 
to 110 mV. These records share most of the features 
seen in Fig. 2 A: namely, delayed activation from the 
pulse onset, increasing rate of rise with membrane de- 
polarizations (Fig. 4 A, inset) , and a lowest voltage at 
which IK V becomes detectable at approximately —20 mV. 
In addition, Fig. 2 A shows that the decay phases of IK V 
records are significandy more pronounced than those 



DiFranco et al. 



115 




-100 5) 0 50 100 -100 -50 0 50 100 

Membrane Potential (mV) Membrane Potential (mV) 



Figure 4. IK V records in osmotically treated fibers and isradipine effects. (A) IK V records elicited by 500-ms depolarizations from Vh to 
membrane potentials spanning from —30 to 90 mV every 20 mV. The inset shows the onset of the currents in an expanded timescale. 
(B) Effects of 1 pM isradipine on IK V . Same fiber and protocol as in A. The onset of the currents is shown in an expanded timescale in the 
inset. (C) Peak IK V (closed circles) and peak gK v (open circles) determined from the data in A (black circles, control conditions) and 
B (red circles, isradipine) . IK V data points are connected by straight lines. The solid black lines are double Boltzmann fits to the control 
gK v data with the following parameters: A = 3.2 mS/cm 2 ; e = 0.21; VI = - 10 mV; Kl = 6 mV; V2 = 31 mV; K2 = 19 mV. The red solid lines 
are double Boltzmann fits to the isradipine gK v data with the following parameters: A =2.4 mS/cm 2 ; e =0.2; VI = —2 mV; Kl = 10 mV; 
V2 = 31 mV; K2 = 19 mV. (D) Peak gK v plots for a population of fibers maintained in control conditions (black circles; 12 experiments, 
1 1 fibers) and a subpopulation of fibers ( n = 4) exposed to 1 uM isradipine (red circles) . Both datasets were fitted with double Boltzmann 
equations. Parameters for control data (black line) were (mean ± SD) A = 2.93 ±0.5 pA/cm 2 ; e = 0.29 ± 0.06; VI = -11.8±3.8mV; Kl = 
5.8 ± 0.7 mV; V2 = 26.5 ± 6.2 mV; K2 = 16.7 ± 3 mV. For this dataset, the fit with a double Boltzmann is >8 x 10 7 times better than with a 
single Boltzmann function. Parameters for isradipine data (red line) were A = 2.17 ± 0.7 mS/cm 2 (*); e = 0.22 ± 0.07; VI = —10.4 ± 2.5 mV; 
Kl = 7 ± 1 mV (*); V2 = 25.6 ± 4.3 mV; K2 = 16 ± 1.7 mV. The * indicates statistical significance with respect to the control (P < 0.05). 
The AIC test indicates that the fit to the isradipine data with a double Boltzmann is only 52-fold better than with a single Boltzmann 
function. Error bars represent SEM. 



observed in intact fibers (particularly at large depolar- 
izations) . Also, decay of IK V in osmotically treated fibers 
is monotonic, without the apparent contaminations by 
ICa-dependent contributions seen in control fibers. 
This is consistent with the almost complete disappear- 
ance of the slow decaying tail currents at the end of the 
pulses. Consistent with our TPLSM observations and 
capacitance measurements, we submit that the results 
in Fig. 4 A can be explained if deeper regions of the TTS 
were effectively disconnected from the rest of the fiber 
by osmotic shock treatment. It is interesting to note 
that although the voltage dependence of the peak IK V 
for the treated fiber (Fig. 4 C, black circles) is compa- 
rable with that obtained in intact fibers (Fig. 1 C), its 



magnitude is smaller. These features are better illus- 
trated in the peak gK v plot in Fig. 4 C (black circles) . 
Moreover, the voltage dependence of peak gK v still sug- 
gests the contribution of two channels to the total cur- 
rents; this is reinforced in Fig. 4 D by showing that only a 
double Boltzmann equation (black trace) with a (gK v ) max 
of 2.77 ± 0.02 mS/cm" (significantly smaller than the 
5 mS/ cm 2 in intact fibers; P < 0.05) can adequately fit 
the results (Fig. 4 C, black open circles) obtained from 
several osmotically treated fibers (n = 7). Because the 
peripheral TTS remnants (together with the sarcolemma) 
are expected to be under more adequate voltage com- 
mands, as compared with the entire TTS in intact fi- 
bers, the data in Fig. 4 C support the notion that two 



116 Potassium currents in mammalian skeletal muscle fibers 



populations of channels (with different voltage depen- 
dence) contribute to the IK V records. 

To determine whether the effects of isradipine de- 
pended on the presence of ICa, we next exposed this 
osmotically treated fiber to 0.5 uM isradipine. Interest- 
ingly, we found (Fig. 4 B) that isradipine still blocks IK V in 
osmotically treated fibers. The effectiveness of the isra- 
dipine block on IKy is quantitatively evaluated at every 
voltage in the peak IK V and peak gK v plots in Fig. 4 C 
(red circles and traces) ; it can be observed that addition 
of isradipine resulted in an overall reduction of the cur- 
rents and conductances at all the voltages tested. Fig. 4 D 
shows that on average, the maximal peak gK v is signifi- 
candy reduced (by 26%) in the presence of 0.5 uM is- 
radipine (black and red circles) . By fitting the peak gK v 
data under control conditions and in the presence of 
isradipine with double Boltzmann equations (Fig. 4 D, 
black and red traces, respectively), we found that (as in 
intact fibers) the blocker tends to affect more promi- 
nently the low threshold component of the currents. 
Nevertheless, probably because of the relatively small 
number of experiments in which the blocker was tested 
(n = 4), the difference in e between control and treated 
fibers was not significant in this case (P = 0.1) ; however, 
the improvement attained by fitting the data to a dou- 
ble Boltzmann instead of a single one is not as marked 
as under control conditions (only 52-fold vs. >10 5 -fold) , 
thus insinuating a lesser contribution of one channel 
component to the total current in isradipine. 

Kinetic analysis of the current components contributing 
to IK V in osmotically treated fibers 

For the reasons discussed in the previous section, cur- 
rent records from partially detubulated fibers are useful 
to obtain (with better accuracy than in intact fibers) the 
voltage dependence of kinetic parameters needed to 
simulate the presence of components contributing to 
IK V . Although we tried initially to predict the properties 
of IK V records using kinetic models contemplating only 
one channel entity, our attempts to reasonably repro- 
duce their major features (delayed activation, onset 
kinetics, inactivation, and overall voltage dependence) 
were unsuccessful. Consequentiy, we focus on the task of 
explaining the properties of IK V in the fiber described 
in Fig. 4, together with data from eight other osmoti- 
cally treated fibers, assuming that two channels (channels 
A and B) contribute to IKy. For simplicity, our models 
(described in the Appendix) assumed that K channel 
activation followed a sequential scheme including four 
closed states (C1-C4) and one open state (O); as it is 
well known, this allows for the delayed activation of 
the currents akin to the n 4 of the Hodgkin and Huxley 
(HH) formulation (Hodgkin and Huxley, 1952b) and 
is compatible with the sequential activation of the 
a subunits in tetrameric K channels. Another simplify- 
ing approximation of our model is that the n 4 -type 



voltage-dependent activation is sequentially coupled 
in a voltage-independent fashion with either a single 
inactivated state (e.g., N type, channel A) or two sequen- 
tial N-type and C-type inactivated states (channel B). 
Examples of our resulting simulations are illustrated 
in Fig. 5 A, and the mean model parameters are in- 
cluded in Table Al of the Appendix. As an example of 
the goodness of the fit with the two-channel model, the 
data in Fig. 5 A (black traces) are superimposed with 
model simulations (red traces) ; it can be seen that ex- 
perimental and simulated currents can be hardly distin- 
guished from each other. The inset in Fig. 5 A illustrates 
that not only the overall magnitude and decay charac- 
teristics of the experimental records, but, most impor- 
tantly, the activation properties of IK V , are fairly well 
reproduced by model predictions. It is interesting to 
note at this point that, to attain this kind of fit to the 
experimental data, it was necessary to include a low 
threshold fully inactivating channel (channel A; Fig. 5 B, 
blue traces) and a high threshold channel that only par- 
tially inactivates (channel B; Fig. 5 B, green traces). In- 
activation was a required feature for both channels. As 
expected, the distinct voltage dependence of activation 
and rate of inactivation between channels A and B de- 
termines that the fraction of the total IKy contributed 
by each channel was voltage and time dependent. For 
example, the inset in Fig. 5 A illustrates that channel A, 
by being responsible for most of the current in response 
to small depolarizations, actually determines the overall 
delay in activation of IK V . In contrast, the rapid activa- 
tion kinetics of channel B dominates the IK V records for 
depolarizing pulses >120 mV. Furthermore, the widely 
different inactivation features of channel B predict that 
current contributions during long depolarizations arise 
exclusively from channel B. The voltage dependence of 
peak IKy and gK v for the total currents in Fig. 5 A are 
shown in Fig. 5 C (black circles) and Fig. 5 D (black 
circles), respectively. The continuous traces are model 
predictions with the parameters listed in Table 2. Like- 
wise, the corresponding fractions of IKy and gK v con- 
tributed by channels A (blue traces) and B (green traces) 
are shown superimposed in the plots of Fig. 5 (C and D) . 
Fig. 5 C shows that, as expected, currents through chan- 
nel A at large depolarizations reach an asymptotic de- 
pendence (not depicted) that intersects the abscissa at 
approximately —90 mV (E K ) . Finally, Fig. 5 D shows that 
channel A reaches a maximal conductance of ^T mS/cnr 
at ^20 mV, whereas the conductance for channel B still 
does not reach a plateau (^2 mS/ cm 2 ) at potentials as 
high as 90 mV. 

Biochemical identification of K v channels in mouse muscle 
In an attempt to determine what molecular entities 
could possibly be responsible for the IKy in adult FDB 
fibers, we performed immunoblotting experiments aimed 
to detect inactivating isoforms of known voltage-dependent 
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Figure 5. Predictions of IK V records from an osmotically treated fiber by a two-channel model. (A) Superposition of experimental 
(black traces) and simulated currents (red traces). Currents were elicited by depolarizations to membrane potentials ranging from 
— 10 to 90 mV in 20-mV increments. The inset displays the onset of currents in an expanded timescale. (B) Individual model current 
contributions attributed to channel A (blue traces) and channel B (green traces) . For every voltage, the currents in A are the sum of the 
traces in B. The inset displays the currents predicted for channels A and B in an expanded timescale. (C) Voltage dependence of peak 
IK V from experimental data (circles) and from currents predicted for channels A (blue trace) and B (green trace). The black trace is 
the sum of the currents from channels A and B. (D) Voltage dependence of peak gK v for the experimental currents (circles) and for 
the currents thought to be A (blue trace) and B channels (green trace) . The black trace was the predicted conductance from the overall 
model simulations. The model parameters for this particular fiber are listed in Table 2 (see Appendix for nomenclature). 



K channels. To this end, MLs and microsomal prepara- 
tions were tested with antibodies raised against shaker 
(K v l. 4) and shaw (K v 3.4) channels. Two different anti- 
bodies recognizing epitopes located at the N terminus 
(75-010; Neuromab) and C terminus (APC-007; Alo- 
mone) of K v 1.4 were used; Fig. 6 A shows that both of 
them recognize one protein band of ^97 kD in total 
protein MLs. Furthermore, both antibodies also report 
one band of comparable molecular weight in MM prepa- 
rations; the slight difference in migration could be ex- 
plained by the enhanced presence of SERCA1, a 1 10-kD 
membrane protein that is highly enriched in skeletal 
muscle MM preparations. The aforementioned results 
are in agreement with similar data reporting that APC- 
007 recognizes an ^97-kD band in heart muscle prepara- 



tions from wild-type animals, which is absent in those 
from K v 1.4-null mice (Guo et al., 2000). Because the 
presence of Kyi. 4 has been extensively demonstrated 
in mouse hippocampus (Sheng etal., 1992, 1993; Maletic- 
Savatic et al., 1995; Rhodes et al., 1997; Monaghan etal., 
2001), as a further control for antibody specificity, we 
also tested lysates from this tissue with the same anti- 
bodies as in muscle preparations. Fig. 6 A shows that 
both 75-010 and APC-007 recognize two bands in the 
HL lanes: a prominent ^95-kD band and a lighter 
band of ^85 kD. This is in agreement with previous 
studies showing the same two bands in mouse hippo- 
campus tissue (Veh etal., 1995; Wickenden et al., 1999; 
Juiz et al., 2000). Altogether, our results and those of 
Guo et al. (2000) may suggest contrasting differences 
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TABLE 2 



Parameter values for K v channels A and B in Fig. 5 
Parameter Channel A Channel B 





0.0232 


0.0155 


A (ms ) 




(ms" 1 ) 


0.01 


0.38 






V * (mV) 


4 


6.9 






Kx (mV) 


14 


16.7 


k BX (mV) 


14.5 


38 


*x (ms -1 ) 


0.004 


0.0074 


<Px (ms -1 ) 


1 x 10 -5 


0.0114 


7b (ms -1 ) 




0.0033 


8 B (ms -1 ) 




1 x 10 -6 


gK^ (mS/cm 2 ) 


1.23 


2.28 


Percent contribution 


0.35 


0.65 



in K v l .4 processing between muscle and brain tissues. 
Western blot analysis also provides supporting evidence 
for the presence of K v 3.4 in FDB muscles. Fig. 6 B 
shows the results obtained when ML and MM prepara- 
tions were tested with a monoclonal anti-K v 3.4 anti- 
body (raised against amino acids 175-192 of this channel; 
75-112, Neuromab). In this case, the antibody detects 
a single band of ^95 kD in the ML lane and, similar to 
what was observed for K v 1.4, a slighdy lighter band in 
the MM preparation. Our finding that K v 3.4 is ex- 
pressed in leg and foot muscles is in agreement with 
a previous study demonstrating the presence of this 
channel in crude membranes of rat sartorius muscles 
(Abbott etal., 2001). 

Simultaneous detection of di-8-ANEPPS transients 
and IKy in intact FDB fibers 

The results from osmotically treated fibers suggested 
that a substantial fraction of IK V arises from the TTS; 
nevertheless, quantitating this fraction is difficult be- 
cause (as shown in Fig. 3 and related text) detubu- 



lation by osmotic shocks is incomplete, leading to 
uncertainties on the extent of TTS contributions that 
still remain in current records. Alternatively, as we 
have extensively demonstrated previously (Heiny et al., 
1983; Ashcroft et al., 1985; DiFranco and Vergara, 
2011; DiFranco et al., 2011a), the relative distribution 
of a given type of ion channel between the surface 
and TTS membrane compartments can be quanti- 
tatively determined by studying the effects that the 
current carried by these channels has on the TTS 
membrane potential, as assessed by potentiometric 
dyes. Fig. 7 shows, in an experiment done in the pres- 
ence of 5 pM nifedipine (which has an effect similar 
to that of 0.5 pM isradipine) to avoid ICa contami- 
nation, that the features of di-8-ANEPPS transients 
are actually altered by the activation of IK V . Fig. 7 A 
demonstrates that small depolarizations (to less than 
-50 mV) do not significantly activate IK V (Fig. 7 C, 
black trace) and result in a step-like fluorescence sig- 
nal (Fig. 7 A, black trace) . This optical response is in- 
distinguishable from that obtained from the same fiber 
after blocking IK V with TEA (Fig. 7 B, black trace). 
Likewise, optical transients recorded in response to 
hyperpolarizing pulses under control or TEA solutions 
could not be distinguished from each other (unpub- 
lished data) . In contrast, depolarizations to 70 mV, for 
instance, which activate a large IK V (Fig. 7 C, orange 
trace), result in di-8-ANEPPS transients that are de- 
pressed early after the pulse onset and later rise in am- 
plitude toward the end of the pulse (Fig. 7 A, orange 
trace), thus creating a notch (or valley) soon after the 
onset of the transient. These features are better noticed 
with an expanded timescale as in the inset of Fig. 7 A 
(same trace colors) . It should be noted that the time 
course of the increase in fluorescence during the pulse 
closely resembles that of the decay phase of the IK V 
trace, clearly suggesting that the activation of the 
latter may be responsible for the kinetic features of 
the di-8-ANEPPS fluorescence transient. That this is 
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Figure 6. Western blot analysis of K v 1.4 and K v 3.4 expression. (A and B) Immunoblots tested with the specific antibodies anti-K v 1.4 (A) 
and anti-K v 3.4 (B). See Materials and methods and Results for details. The results with the 75-010 antibody were obtained from two dif- 
ferent blots, as indicated by the vertical line. The results with the APC-007 antibody were obtained from a single blot, but a lane between 
MM and ML was removed, as indicated by the vertical line. 
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Figure 7. Effects of IK V activation on di-8-ANEPPS transients. (A) di-8-ANEPPS transients recorded in response to 500-ms depolariza- 
tions from Vh to membrane potentials spanning —90-70 mV in 40-mV steps. The inset shows the early phases of fluorescence transients 
in an expanded timescale. (B) Fluorescence transients elicited in the same fiber as in A, and using the same pulse protocol, after ex- 
changing NMG-Tyrode with TEA-Tyrode. (C and D) The currents recorded simultaneously with the optical records in A and B are shown 
in C and D, respectively. For both fluorescence and current records, green, black, red, blue, and orange traces correspond to responses 
to depolarization to —90, —50, —10, 30, and 70 mV, respectively. 



the case is demonstrated by abolishing IK V with TEA. 
Under these passive conditions, the largest depolariza- 
tion elicits a step-like fluorescence response (Fig. 7 B, 
orange trace). It should be noted that the transients 
recorded under control and passive conditions (Fig. 7, 
A and B, orange traces) differ not only in their kinetic 
features, but also in the final value attained at the end 
of the pulse; the control transient reaches a signifi- 
cantly smaller value as compared with that in the pres- 
ence of TEA. These differences can be explained by the 
relatively large IK V still present at the end of the pulse. 
Comparative inspection of Fig. 7 (A and B) clearly shows 
that depolarizations to intermediate values between 



—60 and 110 mV are associated with fluorescence tran- 
sients displaying intermediate characteristics between 
those seen in the black and orange traces. By compar- 
ing traces in Fig. 7 (A and B), it can be noticed that 
the activation of IK V also affects the properties of di-8- 
ANEPPS records after the end of the pulses. For ex- 
ample, the orange and blue traces in Fig. 7 A remain 
substantially elevated above baseline after the end of 
the voltage pulse. It is also clear that these postpulse 
elevations depend on the magnitude of IK V because it 
becomes smaller for smaller pulses (Fig. 7 A, red trace) 
and is effectively eliminated by blocking IK V with TEA 
(Fig. 7 B). An obvious interpretation of these latter 
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Figure 8. Radial cable predictions of IK V and TTS membrane potential changes in intact fibers using multistate models of K channels. 

(A) di-8-ANEPPS transients recorded in the presence of NMG-Tyrode. The fiber was depolarized in 20-mV steps from —90 (V H ) to 110 mV. 

(B) Model predictions of the mean TTS voltage (see Appendix) for the same conditions as in A, and assuming the same value of gK v for 
the sarcolemma and TTS membranes (0.97 mS/cm 2 ) . The model parameters for channels A and B are those in Table Al of the Appen- 
dix. (C) Model predictions of the mean TTS voltage, assuming that gK v (3.88 mS/cm 2 ) is located exclusively at the surface membrane. 
(D) IK V recorded simultaneously with the optical data in A. (E and F) The predicted currents corresponding to the optical simulations 
in B and C are shown in E and F, respectively. The same color code is used for electrical and optical data. 



results is that the activation of IK V results not only in a 
sustained flow of K ions across the membranes of the 
TTS but also leads to a substantial accumulation of 
K ions in the lumen of the T tubules. The concomitant 
changes in E K result in sustained depolarizations of the 
TTS membranes that are accurately reported in poten- 
tiometric dye records for a long time after the com- 
mand pulses are turned off. 

Model simulations of voltage changes in the TTS 
associated with the activation of K v channels 

We have previously shown that when radial cable model 
predictions of TTS voltage changes are contrasted with 
potentiometric dye transients acquired simultaneously 
with the currents in voltage-clamped muscle fibers, 
there are enough constrains in the model to allow for 
quantitative evaluations of conductance distributions 
between the surface and TTS membrane compartments 
(Ashcroft et al., 1985; DiFranco and Vergara, 2011; 
DiFranco et al. , 20 1 1 a) . The specific model used here to 
simulate IK V and di-8-ANEPPS is described in detail 
in the Appendix. In brief, it is similar to that used pre- 
viously for the cases of C1C-1 and Nayl -4 channels but 
includes instead the multistate models of two K channels 



(refined with data from osmotically treated fibers) and 
equations to account for K accumulation. In addition, 
the relative contribution of radial shells of TTS mem- 
brane to the experimental fluorescence transients takes 
into account a correction that incorporates the optical 
properties of the 100x, 1.4 numeral aperture objective 
used, as previously described (Kim and Vergara, 1998) . 
The parameter values used in the model are given in 
Tables Al and A2 of the Appendix. Fig. 8 (A and D) 
shows optical transients and IK V recorded in the ab- 
sence of Ca blockers, respectively, in response to a fam- 
ily of 25-ms duration depolarization pulses. It can be 
seen in Fig. 8 A that the optical transients display the 
typical current-dependent attenuations during the pulse, 
as described in Fig. 7 A. It can also be observed that at 
the end of the pulse, the fluorescence signals decay 
slowly to baseline, in correlation with the relatively slow 
decaying tail current observed in Fig. 8 D. For the same 
depolarizations, the simulated mean TTS membrane 
potential (corrected for optical detection; see Appendix) 
is shown in Fig. 8 B; the simultaneous model predic- 
tions of IKy are shown in Fig. 8 E. It can be seen that the 
simulated traces quantitatively reproduce most of the 
kinetic features of both sets of experimental data during 
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each of the depolarizations tested. Nonetheless, the 
model does not reproduce the slow decay of the optical 
transients and the tail currents seen experimentally 
after the end of the pulses; instead, model data show 
smaller tail currents and postpulse depolarizations. 
These discrepancies can be explained if we assume 
that optical records and tail currents result from re- 
sidual IK 1R (not sufficiently blocked by 5 mM Rb) and/ 
or from unblocked Ca 2+ currents that are not included 
in the model. Also, the model does not contemplate the 
existence of charge movement currents present in the 
experimental records. It is important to note that, al- 
though the model predicts substantive K accumulation 
in the lumen of the TTS occurring concomitantly with 
IK V during the pulses, these are not manifested as tail 
currents and/or postpulse depolarizations. Another no- 
ticeable difference between model predictions and data 
is that for large depolarizations, theoretical IK V traces 
decay slighdy faster than the experimental ones. 

We may now ask the hypothetical question of what the 
predicted TTS voltage changes and currents would look 
like had all gKy been located at the sarcolemma. To this 
end, as shown in Fig. 8 F, we account reasonably for the 
total IK V by including a large gK v (4 mS/ cm 2 ) at the 
sarcolemma. However, as shown in Fig. 8 C, the mean 
TTS membrane potential changes display step-like re- 
sponses for all depolarizations that are in sharp contrast 
with the observed optical responses. Notice that the pre- 
dicted TTS transients do not rise instantaneously as a 
result of the relatively slow charging of the large capaci- 
tance of the TTS, but their amplitude depends linearly 
on the amplitude of the command pulse. It should be 
also observed that the responses shown in Fig. 8 C are 
quite similar to those obtained in the presence of TEA 
(e.g., Fig. 7 B). Finally, as expected, no tail currents are 
predicted by the model when currents do not flow 
across TTS membranes. 

Quantitative evaluation of current-dependent attenuations 
on di-8-ANEPPS transients 

The effects of IK V on TTS membrane potential changes, 
as assessed by di-8-ANEPPS transients, are further ex- 
plored in Fig. 9, which shows superimposed optical 
records obtained in response to two pulse amplitudes 
(100 and 180 mV; Fig. 9, A and B, respectively) before 
(red traces) and after (black traces) blocking the cor- 
responding IKy, as shown in Fig. 7 (C and D; same color 
coding) , respectively. The IK V associated with the 100-mV 
depolarization generates an observable reduction in 
the magnitude of the optical trace (Fig. 9 A, red trace), 
which is larger early during the pulse than at the end of 
the pulse. For the larger pulse, the attenuating effects of 
the larger current (Fig. 9 D, red trace) on the optical 
record (Fig. 9 B, red trace) are more prominent early 
during the pulse, but they remained throughout the du- 
ration of the pulse. It should be noted that the magni- 



tude of the depression displays a reversed pattern with 
respect to the decay of the corresponding IK V record. 
To quantitatively evaluate the relationship between IK V 
activation and the transient depression that is observed 
in di-8-ANEPPS transients (creating a valley or notch) 
relative to the step-like shapes seen in the absence of 
current with external TEA, we must implement a nor- 
malization procedure. Furthermore, this normalization 
is necessary to later compare the optical data with model 
predictions that compute the mean membrane potential 
changes in the TTS. To this end, we define attenuation 
as the percent difference between transients recorded 
before (Fig. 9, A and B, red traces) and after blocking 
IK V (Fig. 9, A and B, black traces), normalized by the 
amplitude of the largest transient acquired (e.g., the black 
trace in Fig. 9 B for this particular experiment). As in 
this experiment, the normalization is usually done with 
respect to signals elicited (after making the fibers elec- 
trically passive) by 200-mV pulses, a typical maximum 
stimulus. Attenuation traces calculated this way are shown 
in green in Fig. 9 (A and B; and in the inset of Fig. 9 B) . 
Several features are salient: (a) their amplitude is sig- 
nificantly larger for the larger depolarization (~9% for 
Fig. 9 A and ^26% for Fig. 9 B) , in correlation with the 
magnitude of the IK V records shown in Fig. 9 (C and D) ; 
(b) the time course of the attenuation traces mirrors 
that of the corresponding IKy records; namely, it has 
rising and decay phases comparable with those of the 
current traces; and (c) the attenuation trace reaches a 
peak with a visible delay with respect to the onset of the 
pulse, and its peak value (peak attenuation) closely co- 
incides in time (within a millisecond) with the notch 
(minimum) observed in di-8-ANEPPS transients and 
with the peak of IKy. The close correlation between the 
peak of the attenuation trace and the notch in the opti- 
cal transients itself is better visualized at the expanded 
timescale of the inset to Fig. 9 B. Altogether, these fea- 
tures suggest a kinetic correlation between the atten- 
uation of optical records and IKy that could possibly 
be exploited to obtain a mechanistic understanding of 
their interdependence (see the following paragraph) . 

Another manifestation of the current dependence of 
optical traces is illustrated in Fig. 10. The red circles in 
Fig. 10 A correspond to AF/F values of the notches in 
di-8-ANEPPS transients plotted as a function of mem- 
brane potential, whereas the black circles show the 
steady-state value of the transients in the presence of 
TEA. The depression in magnitude of the transients 
when the currents are active is seen as markedly sublin- 
ear voltage dependence; this dramatically contrasts with 
the quasilinear behavior in their magnitude when the 
currents are blocked. For comparison, the voltage de- 
pendence of the peak IKy is shown in Fig. 10 B. It can be 
seen that the departure from the quasilinearity of the 
optical transients in the absence of currents (Fig. 10 A, 
compare black with red circles) occurs precisely in the 



122 Potassium currents in mammalian skeletal muscle fibers 



A 




800 



600 



£ 400 
o 



200 



-200 



100 200 300 400 500 600 700 



200 400 
Time (ms) 



600 



B 




D 



0 100 200 300 400 500 600 700 



800 
600 
400 
200 
0- 
-200 




200 400 
Time (ms) 



600 



Figure 9. Calculation of Independent attenuation of fluorescence transients. The attenuations in di-8-ANEPPS transients are exempli- 
fied with two pulse amplitudes (100 mV in A and 200 mV in B) . (A) The red and black traces are the di-8-ANEPPS transients recorded 
in the presence of NMG-Tyrode and TEA-Tyrode, respectively. The green trace represents the attenuation (percentage), calculated by 
subtracting the red from the black traces and dividing the result by the amplitude of the black trace in B. (B) Same as in A, but for a 
depolarization to 110 mV. The inset shows the initial part of the records in an expended timescale. (C and D) The currents recorded 
simultaneously with the optical data in A and B are shown (with the same color code) in C and D, respectively. To reduce the increment 
of noise inherent to the subtraction operation, the attenuation traces (green traces) were digitally filtered. 



same voltage range at which IK V becomes activated 
(Fig. 10 B, compare black with red circles). This correla- 
tion is reinforced by the fact that the voltage dependence 
of the peak attenuation (Fig. 10 C) closely resembles that 
of the current in Fig. 10 B. Fig. 10 D illustrates compara- 
tively (as in Fig. 10 A) the voltage dependence of notch 
values in di-8-ANEPPS transients before and after IK V 
block obtained from a population of 10 fibers (without 
using Ca channel blockers) ; these conditions allow us to 
obtain an assessment of the role of IK V with physiolo- 
gical values of K conductance. It must be noted that 
the activation of ICa, mainly in response to large depo- 
larizations, does not have significant effects on these 
measurements because the notches in di-8-ANEPPS 
transients occur too early (within <8 ms) to be affected 
by currents that typically peak at >50 ms for the largest 
pulses used here (DiFranco et al., 2011b). Because of 
the larger peak IK V (Fig. 10 B), reaching mean values 



up to 960 uA/cm 2 for depolarizations to 110 mV, the de- 
viation of the notches (from the more linear behavior 
when the currents are blocked) is significantly more promi- 
nent. In this case, the mean AF/F at 110 mVis —0.16, sig- 
nificantly smaller than the —0.24 recorded in TEA The 
close correspondence between the activation voltages for 
IKy and the departure from linearity of di-8-ANEPPS tran- 
sients are clearly indicated by comparing the plots in 
Fig. 10 (D and E). However, this notion becomes further 
stressed by noticing the remarkable similarity between 
Fig. 10 F, which plots the voltage dependence of the peak 
attenuation, and the peak IK V plot in Fig. 10 E. The impor- 
tance of the representation in the form of percent peak 
attenuation is that, as previously demonstrated (DiFranco 
and Vergara, 2011; DiFranco et al., 2011a) and discussed 
later here, it allows for the comparison between model 
predictions of TTS voltage changes (in millivolts) with ex- 
perimental data (in Afarads/farads) . 
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Figure 10. Voltage dependence of the peak attenuation. (A) Red circles: voltage dependence of the minimum (notch) AF/F value 
observed in di-8-ANEPPS transients from a fiber exposed to NMG-Tyrode containing 0.5 uM isradipine. Black circles: steady-state ampli- 
tude of di-8-ANEPPS transients recorded from the same fiber after replacing NMG-Tyrode with TEA-Tyrode. (B) Voltage dependence of 
peak IK V (red circles) and leak current (black circles) measured from the same fiber and conditions as in A. (C) Voltage dependence of 
the percent peak attenuation calculated from data in A. (D) Red circles: voltage dependence of the mean notch AF/F calculated from a 
family of di-8-ANEPPS transients in 10 fibers exposed to NMG-Tyrode (with no added Ca blocker). Black circles: steady-state amplitude 
of di-8-ANEPPS transients recorded from the same population of fibers after replacing NMG-Tyrode with TEA-Tyrode. (E) Voltage de- 
pendence of the mean peak IK V (circles) determined from the same fiber and conditions as in D. (F) Voltage dependence of the mean 
percent peak attenuation calculated from data in D. Error bars represent SEM. 



Quantitative estimation of K v channel distribution 
between the surface and TTS membranes 

The question that has yet to be addressed is what are 
the respective contributions of the sarcolemma and the 
TTS membrane K v currents to the total IK V records 
that are compatible with the optical data? Although the 
question is complex, we already ruled out the possibility 
that all Kv channels are in the sarcolemma because this 
option is generally incompatible with all the optical 
data presented throughout the paper and is specifically 
addressed in Fig. 10 (C and F). The other extreme 
case, that all the K v channels are in the TTS, seems un- 
realistic because it is unlikely that the sarcolemma is 
entirely deprived of a rapid repolarization mechanism; 
also, model simulations suggest that K accumulation in 
the TTS lumen would be excessive if this was the case 
(unpublished data). It might seem that the excellent 
predictions by model simulations of the overall prop- 
erties of optical and current records in the case of 
the fiber in Fig. 10 have already provided a definite an- 
swer about the issue of gK v distribution. However, this 
is not necessarily the case because we should explore 



other possibilities to encompass possible fiber-to-fiber 
variability in the data. Fortunately, our simultaneous 
measurements of the voltage dependence of peak atten- 
uation and peak I-V curves from multiple fibers (Fig. 10, 
E and F) allow us to complete the task. Fig. 11 shows 
the results from model simulations using various (sur- 
face/TTS) gK v ratios to generate a range of attenuation 
profiles while constraining the total peak IK V depen- 
dence to be constant. The resulting simulated I-V plots 
(Fig. 11, D-F, stars) are superimposed with the mean 
experimental data (Fig. 11, D-F, circles) for comparison. 
Fig. 1 1 A shows the superposition of the experimental 
peak attenuation (black circles) with that predicted 
by model simulation assuming a 50:50 (gKy-S/gKy-TTS) 
ratio (red circles). It can be seen that this assumption 
provides a very adequate prediction of the maximum 
attenuation (^30%; Fig. 11 A) and of the maximal 
peak IK V (^960 uA/cm 2 ; Fig. 11 D) obtained from the 
population of fibers. Notably, because the TTS repre- 
sents the largest of the two membrane compartments 
in the muscle fiber, an equal distribution of gK v between 
the TTS and surface membranes results in a predicted 
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Figure 1 1 . Model predictions of overshoot and currents for various gK v distributions between the surface and TTS membranes. 
(A-C) Voltage-dependent peak attenuations calculated from model predictions of the TTS membrane potential using the following 
gKy-S/gKy-TTS ratios: 50:50 (A, red circles), 70:30 (B, green circles), and 30:70 (C, blue circles). The values of gK v were 1.26 mS/cm 2 
for both surface and TTS membranes in A; 2.03 mS/cm 2 at the surface membrane and 0.87 mS/cm 2 at the TTS membranes in B; and 
0.71 mS/cm 2 at the surface membrane and 1.66 mS/cm 2 at the TTS membranes in C. In every panel, the black circles are the mean 
experimental attenuation values shown in Fig. 10 F. (D-F) The mean peak I-V plot shown in Fig. 10 E (black circles) is superimposed 
with model peak current predictions for 50:50 (A, red circles) , 70:30 (B, green circles) , and 30:70 (C, blue circles) gKy-S/ gKy-TTS ratios, 
respectively. In every panel, stars represent the total current, closed triangles represent the TTS current component, and open triangles 
represent the surface membrane current component. The fiber radius used in model simulations was the mean for all the experiments 
(24 pm), and R s was 40 flcm 2 . The rest of the model's parameters are listed in Tables Al and A2 of the Appendix. Error bars rep- 
resent SEM. 



maximal IK V arising from the TTS of ^690 pA/ cm 2 
(at 110 mV) . This represents ^72% of the total current 
at this membrane potential, whereas the surface mem- 
brane contributes the other ^28% (^270 pA/cm 2 ). 
Fig. 11 (B and E) shows the results of simulations in 
which we limited the density of K v channels in the TTS 
membranes to be ^43% of that at the surface membrane 
(70:30 ratio). It can be seen that this ratio predicts that 
the maximal TTS voltage attenuation (Fig. 11 B, green 
circles) is ^23%, significantly smaller than the one 
observed experimentally. Thus, although the peak I-V 
plot of the total current (Fig. HE, stars) was similar 
to the one shown in Fig. 11 D, the TTS component 
(Fig. HE, closed triangles) with a maximum of ^541 
pA/ cm 2 seems to be too small to explain the observed 
peak attenuation in the mean population of fibers. In 
spite of this, this current component would still contrib- 
ute a substantial proportion (>50%) to the total current 
recorded at the same membrane potential. These types 
of comparisons illustrate the great sensitivity afforded 
by the optical measurements to discriminate between 



channel distributions that would give the same total 
current but have very different balances between the 
surface and TTS current contributions. Without the in- 
formation provided by the potentiometric dye records 
about the associated changes in TTS membrane poten- 
tial, it would be impossible to decide which one is correct. 
Likewise, a 30:70 gKy-S/gKy-TTS ratio (Fig. 11, C and F) 
predicts a peak overshoot of ^34% (Fig. 11 C, blue 
circles) , which is larger than those recorded. In this case, 
the current component arising from the TTS (Fig. 1 1 F , 
closed triangles) has a maximum of ^808 pA/cm 2 , which 
not only seems to be too large to explain the optical 
data but would also account for >84% of the total cur- 
rent at that potential. This is also incompatible with 
detubulation experiments. The suggestion from this 
analysis comparing experimental data and model sim- 
ulations at a fixed total current is that the overall distri- 
bution of K v channels between the surface and TTS 
of mammalian skeletal muscle fibers is probably con- 
strained within the range between 40:60 and 60:40 
(gKv-S/gKy-TTS). 
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DISCUSSION 

Though it is well accepted that IK V plays a crucial role in 
muscle excitability by being the principal component 
responsible for the repolarization phase of the AP, there 
is a great shortage of information in the literature about 
the detailed properties of gK v in mammalian skeletal 
muscle fibers. The paucity is even more striking for 
mouse muscle fibers, which is unfortunate because this 
is the animal choice for models of human diseases, in- 
cluding mutations that affect K v channels. Furthermore, 
up until now, there was uncertainty about the presence 
of gK v in the TTS, which is a very relevant issue regard- 
ing the possibility that K accumulation in the TTS 
lumen during bouts of sustained activity would possibly 
lead to depolarization of the muscle fibers, eventually 
compromising their excitability. This paper, by provid- 
ing a quantitative description of the properties of gK v , 
including the involvement of TTS, not only fills in gaps 
of understanding about the physiological properties 
of this conductive pathway in murine skeletal muscle 
fibers but also opens new avenues for the investigation 
of its alterations in diseased conditions. Experience ac- 
cumulated during the last decade suggests that short 
muscle fibers isolated from the toes of adult mice are 
practical and reliable preparations for electrophysiolog- 
ical studies, using a two-microelectrode voltage-clamp 
technique (Woods et al., 2004, 2005; Ursu et al., 2005; 
Lueck et al., 2010; DiFranco et al., 2011a,b; Fu et al., 
2011). We have recently demonstrated that this classi- 
cal electrophysiological approach can be combined 
with optical measurements of TTS membrane potential 
changes and model simulations to determine, with min- 
imal invasion, the quantitative allocation of particular con- 
ductances in the TTS and surface membrane (DiFranco 
and Vergara, 2011; DiFranco etal., 2011a). In this work, 
we use this approach to investigate overall properties of 
IK V in muscle. 

IK V in intact FDB muscle fibers 

The voltage- and time-dependent IK V recorded from 
intact and detubulated fibers display all the canonical 
features that traditionally have been ascribed to K v 
channels: a distinct activation potential, outward rec- 
tification, voltage-dependent delayed onset, voltage- 
dependent rate of rise, and incomplete inactivation. The 
amplitude of IK V recorded from intact fibers was large; 
for instance, for depolarizations similar to those reached 
at the peak of the AP (^45 mV) , peak IK V can reach val- 
ues of ^500 uA/cm 2 , which is smaller, but within range, 
of the maximal Na current for this preparation (unpub- 
lished data) . No quantitative comparative measurements 
of IKy have been provided for mouse muscle fibers be- 
cause most previous work was performed using either 
on-cell or excised-patch configuration of patch clamp 
(Brinkmeier etal., 1991; Hocherman and Bezanilla, 1996) 



or whole-cell patch clamp; the latter method yielded 
K currents estimated to be ^10% of the magnitude of 
those reported in this paper (Brinkmeier et al., 1991). 
In contrast, our mean peak IK V values of up to 990 uA/ cm 2 
at 110 mV are in reasonable agreement with results 
reported for rat muscle fibers using other voltage- 
clamp approaches (Duval and Leoty, 1980; Pappone, 
1980; Beam and Donaldson, 1983a). Also, the main 
kinetic features of our IK V records are similar to those 
recorded from fast muscle fibers of mouse and rat, 
showing a single peak (Duval and Leoty, 1980) and in- 
complete inactivation during 500-ms pulses (Duval and 
Leoty, 1980; Beam and Donaldson, 1983a; Brinkmeier 
etal., 1991; Hocherman and Bezanilla, 1996). A relatively 
large fraction of IK V does not inactivate even for 1-s 
pulses (unpublished data) , but we have observed that 
the inactivation process is cumulative. Furthermore, we 
determined that a resting period of 7-10 s was sufficient 
to avoid this phenomenon; consequendy, the prepara- 
tion was allowed to rest for at least this long between 
pulses during the acquisition of a family of records. 
More importantly, the decay kinetics of IK V records is a 
complex process involving at least two time constants; 
although this may be suggestive of more than one chan- 
nel contribution, we preferred to use a more direct ap- 
proach to investigate this issue. Interestingly, a two-time 
constant decay was also found in ensemble currents of 
membrane patches from mouse fibers (Hocherman and 
Bezanilla, 1996). 

A complication found in our experiments was the con- 
tamination by ICa in IK V records from intact fibers and 
its eventual elimination by the use of typical Ca 2+ channel 
blockers. In reality, ICa's are small relative to the magni- 
tude of IK V , but they do manifest as small downward de- 
flections during the decay phases of IKy records elicited 
by depolarizations larger than 140 mV and longer than 
100 ms (Fig. 2 A) . Prominent ICa tails also dominate the 
current records after the end of the pulses. It could have 
been argued that removal of Ca 2+ ions from the extracel- 
lular solution would have eliminated ICa contamina- 
tions. However, we avoided this approach mainly for two 
reasons: the cells became leaky, probably because of de- 
terioration of the membrane seal around the microelec- 
trode tips, and the voltage dependence of K channels was 
expected to be substantially shifted as a result of changes 
in surface charge screening (Hille et al., 1975; DiFranco 
et al., 201 lb) . We have found that Mg 2+ is a poor replace- 
ment of Ca 2+ toward these ends. A reasonable replace- 
ment for Ca 2+ is Ba 2+ , but it could not be used because it 
blocked IKy and is permeable through Ca 2+ channels. 
When we tried large concentrations of external Cd 2+ 
(>1 mM), we found that the leak currents increased to 
levels incompatible with the experiments (unpublished 
data); also, this manipulation resulted in a significant 
rightward shift of the voltage dependence of IKy. For these 
reasons, we took the approach of testing dihydropyridine 
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derivatives (e.g., nifedipine and isradipine) and other 
organic blockers of Cayl.l. Unfortunately, we found 
that, as extensively described in Results, these drugs 
also block IK V . 

By analogy with other inactivating currents, we stud- 
ied the voltage dependence of IK V by plotting the peak 
currents in response to every pulse as a function of 
the membrane potential reached. We found that peak 
IK V plots follow the typical quasilinear dependence 
for membrane potentials more positive than —10 mV 
(Fig. ID). Assuming a linear driving force (V-E K ) , we cal- 
culated peak gK v and studied its voltage dependence. 
A novel observation is that the peak gK v versus mem- 
brane voltage graphs cannot be accurately predicted by a 
single Boltzmann equation; there is a marked inflection 
at ^25 mV. Instead, the data could be well fitted by dou- 
ble Boltzmann equations with distinct parameters (Fig. 1 , 
B and D; and Fig. 2 D) . The suggestion from these fits is 
that IKy currents contain the contributions of two K v 
channels, one (channel A) with low threshold for activa- 
tion and steep voltage dependence and the second with 
a high activation threshold and less steep voltage de- 
pendence (channel B). The double Boltzmann fit of 
peak gK v plots from intact fibers under control condi- 
tions (Fig. 1 D) suggests a channel A/channel B ampli- 
tude ratio of ^0.4 at the time the current peaks. As 
mentioned previously, peak gK v calculated from peak 
currents are not significantly distorted by ICa con- 
taminations because this current activates very slowly 
(DiFranco et al., 2011b); for example, at 4 ms, when 
IK V peaks for a pulse to 1 10 mV (Fig. 1 A) , we have mea- 
sured that ICa is not larger than 8 uA/cm (not depicted) , 
which is negligible compared with the >900-pA/ cm 2 
peak IKy. It must be noted that, although the data are 
well fitted with double Boltzmann equations, we can- 
not discard the possibility that more than two channels 
actually contribute to IK V . Nevertheless, our findings 
in fibers isolated from the FDB muscle, a typical fast 
type in the mouse (Gonzalez et al., 2000), differ from 
reports that in fast rat fibers, K currents are contrib- 
uted by one channel type (Duval and Leoty, 1980; Beam 
and Donaldson, 1983a); instead, they are compatible 
with the case in soleus (slow) fibers that suggest two com- 
ponents (Duval and Leoty, 1980). 

IK V from osmotically treated fibers 

Our laboratory and another group have recently shown 
that the osmotic shock treatment, by disconnecting a 
significant fraction of the TTS from the surface mem- 
brane, can be used to generate a simpler muscle model 
preparation to identify the properties of the Na con- 
ductance under voltage-clamp conditions (DiFranco 
and Vergara, 2011; Fu et al., 2011). The results pre- 
sented here strongly reinforce this concept, but for K v 
channels. It must be kept in mind, however, that in our 
hands, the formamide osmotic shock procedure yields 



fibers from within the same batch treatment with variable 
degrees of detubulation; furthermore, even within the 
same fiber, there are zones showing more TTS stumps 
than others (unpublished data). We also found that, 
judging by the capacitance and di-8-ANEPPS staining, 
detubulation was never complete. For these reasons, we 
routinely ascertained that every electrophysiological 
experiment was performed in fibers affording the best 
degree of TTS disconnection. For this, we carefully se- 
lected fibers showing di-8-ANEPPS staining only at the 
periphery throughout most of their length (Fig. 3) 
and with capacitances <2.6 pF/cm 2 . In the best cases, 
peripheral rings of TTS remnants spanning 2-3 pm in 
thickness were still observed, a result that explains why 
even this selected population of fibers had a mean ca- 
pacitance of 2.25 pF/ cm 2 , which is significandy larger 
than the theoretical ideal of 0.9 pF/crrr ascribed to 
the sarcolemma (see Appendix) . 

A clear effect of the osmotic shock treatment is a sig- 
nificant reduction of ^40% (P < 0.05) in the peak IKy 
and gK v values compared with those in intact fibers, 
with no apparent impairment of other electrical fea- 
tures. This result unquestionably indicates that a fraction 
of K v channels is located in regions of the TTS discon- 
nected by the procedure. What it does not necessarily 
imply is that only 40% of the current arises from the 
TTS compartment; this will be discussed in the section 
Quantitative evaluation of the relationship between 
attenuation and currents. Another difference between 
IKy recorded from control and treated fibers is the sig- 
nificant reduction in the tail currents in the latter, which 
is expected because Cayl.l channels are mostly located 
in the TTS. Aside from these differences, we found that 
IKy from osmotically shocked fibers (Figs. 4 and 5) dis- 
play features similar to those from intact fibers (Figs. 1 
and 2), which readily indicates that at least a significant 
portion of the decay in the IKy records is a bona fide in- 
activation process; namely, IKy in fibers with a signifi- 
cantly reduced diffusional volume (TTS lumen) , and 
thus with a reduced possibility for K accumulation, still 
display substantial decays. These results are concurrent 
with those obtained using the on-cell configuration of 
patch clamp (Brinkmeier et al., 1991) in which currents 
recorded from surface patches of membrane still show 
inactivation with similar kinetics as those shown here. 
In addition, the properties of IK V in osmotically treated 
fibers seem to be more compatible with those described 
for fibers from the soleus than the iliacus (fast) muscles 
in the rat (Duval and Leoty, 1980). 

A two-channel model of IK V in FDB fibers 

The voltage dependence of peak gK v from detubulated 
fibers also provided valuable information. First, as was 
the case in control fibers, double Boltzmann equations 
were required to appropriately fit the data. This is an 
interesting result because it demonstrates that the two 
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putative channels proposed to explain the data in in- 
tact fibers are not exclusively located in the TTS mem- 
branes, but instead they may be ubiquitously distrib- 
uted. In fact, we took advantage of the results in these 
preparations to obtain the typical behavior of K v chan- 
nels in an ideal single compartment model. Although 
this is only an approximation, we reasoned (for reasons 
already discussed) that gK v in osmotically treated fibers 
arise from the sarcolemma and peripheral regions of 
the TTS, where the membrane potential is much closer 
to that controlled by the voltage-clamp system. 

Altogether, the two-channel model described in 
the Appendix, by quantitatively predicting the voltage- 
dependent kinetic features of data from treated fibers 
(Fig. 5), allows us to gain confidence in the proposal 
that, in skeletal muscle fibers, IK V must emerge from at 
least two channel contributions, one of which does not 
inactivate completely. To predict this behavior, we used a 
general sequential model with the traditional four closed 
and one open state followed by open channel inactiva- 
tion processes as described elsewhere (DeCoursey, 1990; 
Demo and Yellen, 1991; Hoshi et al., 1991; Bett et al., 
2011). Although it was conceivable that both channels 
displayed incomplete (N and/or C) inactivation, the 
adjustment of the model to the data suggested (in all 
fibers analyzed) that the scheme illustrated in Fig. 5 
(channel A: inactivating, low threshold, and steep volt- 
age dependence; and channel B: incompletely inacti- 
vating) was the simplest and strongest. As expected, the 
voltage dependence of peak gk v in channels A and B 
provides a mechanistic explanation for the detailed 
features of peak gK v plots (requiring double Boltzmann 
fits) both in osmotically treated and intact fibers. 

Experiments in osmotically treated fibers also demon- 
strated that, as in intact fibers, IK V is blocked by the or- 
ganic Ca channel blockers isradipine and nifedipine 
(unpublished data) . The data suggest that both putative 
channels (channels A and B) are sensitive to isradi- 
pine because double Boltzmann fits were required when 
fitting peak gK v data in intact and osmotically treated 
fibers. Nevertheless, as shown by the data in Fig. 2 D 
and suggested in Fig. 4 D, channel A seems to have 
an increased sensitivity to the drug. Although further 
studies will be required to characterize the actual block- 
age of IKy by these agents and their specific effects on 
every channel's contribution to IK V records, our results 
from both intact and osmotically treated fibers bear 
important warnings regarding the use of dihydro- 
pyridines when studying the properties of K v channels 
in skeletal muscle and for potential side effects con- 
comitant to their routine use to treat patients with car- 
diac problems. 

The value of an empirical model capable of predict- 
ing the kinetic properties of IK V is double; as stated 
previously, it provides a mechanistic view of what could 
be reasonable contributions of individual channels to 



IK V records. In addition, from a practical point of view, 
it allows us to study the potential implications that the 
presence of gK v has on the membrane potential of the 
TTS when the K channel model equations are incor- 
porated into the radial cable model as described in 
the Appendix. 

Western blotting 

With the idea in mind that at least two functional K v 
channels may be contributing to IK V records, we thought 
it necessary to complement the electrophysiological 
data with biochemical evidence on the actual expres- 
sion of K v channels in FDB and interosseous muscles. 
Although our results are not exhaustive (more channels 
could be expressed) , the Western blot results presented 
in Fig. 6 convincingly show that the muscles used in 
electrophysiological studies express both K v 1.4 and 
K v 3.4. In fact, we also preliminarily tested the expres- 
sion Kyi. 5 and found it in crude extracts (unpublished 
data) . We specifically focused on the expression of 
Kyi. 4 and K v 3.4 because it is generally suggested that 
they might be abundant in skeletal muscle; neverthe- 
less, our results seem to be the first proof that they are 
actually expressed in significant quantities in muscles 
from the mouse. Several lines of evidence support this 
claim regarding Kyl.4: (a) two different antibodies raised 
against epitopes at the N terminus (75-010) and the 
C terminus (APC-007) of this channel identified a unique 
band in both crude extracts and microsomal prepara- 
tions from muscle; (b) the specificity of one of these 
antibodies (APC-007) was previously demonstrated in the 
heart by comparing results from wild-type and K v 1.4-null 
mice (Guo etal., 2000); (c) our results show that APC-007 
recognizes the same — 97-kD band in skeletal muscle as 
that reported for wild-type heart muscle (Guo et al., 
2000); (d) the 75-010 antibody recognizes an epitope 
at the N terminus of K v 1.4 that is exclusive to this pro- 
tein, significantly reducing the chance of cross-reactivity 
with a subunits of the K v l family (Juiz et al., 2000) ; this 
was demonstrated by Neuromab); and (e) the — 97-kD 
band detected by both antibodies in muscle prepara- 
tions coincides with one of the two bands reported in 
mouse hippocampus extracts, a tissue in which the ex- 
pression of K v l. 4 is generally established (Sheng et al., 
1992, 1993; Maletic-Savatic et al., 1995; Veh et al., 1995; 
PJiodes et al., 1997; Wickenden et al., 1999; Juiz et al., 
2000; Monaghan et al., 2001). The —97 kD reported in 
this paper for K v l .4 is not only consistent with that esti- 
mated from Western blots in other tissues but is also 
larger than that stricdy deduced from its amino acid se- 
quence (74 kD); this discrepancy is likely caused by in 
vivo glycosylation of the protein (Sheng et al., 1992; Veh 
et al., 1995; Juiz etal., 2000). 

The evidence for the expression of K v 3.4 in skeletal 
muscle is not as persuasive as for Kyi. 4 but is sustained 
mostly by the specificity of the monoclonal antibody 
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75-112, which has been also demonstrated not to 
cross react with other K v channels (Neuromab). As an 
additional evidence of the specificity of this antibody, 
it has been used to detect changes in the expression 
of K v 3.4 in hippocampus of epileptic rats (Pacheco 
Otalora et al., 2011). Interestingly, in agreement with 
our results, a band of ^95 kD has also been identified 
(by another antibody) in rat sartorius muscle prepara- 
tions (Abbott et al., 2001). As with K v 1.4, the esti- 
mated molecular mass of K v 3.4 is larger than the 70 kD 
predicted from its amino acid sequence (Schroter et al., 
1991); the difference has been explained previously 
on the basis of extensive glycosylation (Cartwright 
etal., 2007). 

Although the biochemical evidence is solid, we do not 
necessarily imply that there is a strict correspondence 
between the properties of the two channels (A and B) 
predicted from our electrophysiological data and those 
of the channels identified by Western blotting. Against 
this possibility, it has been shown that when Kyi. 4 and 
K v 3.4 are expressed in heterologous systems, the re- 
corded currents display properties not identical to 
either channel A or B (Ruppersberg et al., 1990; Po 
et al., 1993; Abbott et al., 2001). However, known fac- 
tors may contribute to the discrepancy: (a) the expres- 
sion of channels in heterologous systems may result 
in functional properties widely different from those 
seen in the muscle fibers, as already shown for the 
Nayl-4 channel (DiFranco and Vergara, 2011; Fu et al., 
201 1) ; and (b) it has been shown that K v channels form 
heterotetramers (Ruppersberg et al., 1990; Po et al., 
1993), associate with regulatory subunits (Abbott et al., 
2001; Pongs and Schwarz, 2010), and are under tight 
regulation in vivo (Covarrubias et al., 1994). As a con- 
sequence, it is difficult to predict the voltage depen- 
dence, kinetics, and pharmacology of native K v channels 
in fully developed muscle fibers from data obtained in 
expression systems (Ruppersberg et al., 1990; Po et al., 
1993; Hashimoto et al., 2000; Abbott et al., 2001). In 
fact, that this may be the case is suggested from expres- 
sion experiments combining more than one K v isoform 
and/or regulatory subunit (Lee et al., 1996; Grunnet 
etal., 2003). Overall, these limitations emphasize that 
the full characterization of IK V while preserving the 
normal cellular context, e.g., under the structural and 
regulatory conditions as they happen in normal cells, 
represents a valuable approach. 

di-8-ANEPPS transients' association with IK V 

As shown in Fig. 7, depolarizations able to activate 
prominent IK V are associated with the appearance of 
prominent depressions in di-8-ANEPPS transients. The 
overall reduction in magnitude of optical signals in 
the presence of IK V with respect to those when the cur- 
rents are blocked by TEA illustrates that activation of 
gK v in the TTS membranes leads to a mean decrease 



of the TTS membrane potential changes with respect 
to those imposed by the voltage clamp at the surface 
membrane. Similar effects were reported previously 
for amphibian muscle fibers as a result of the activa- 
tion of K IR currents (Heiny et al., 1983; Ashcroft et al., 
1985) and for mammalian fibers in association with 
C1C-1 currents (DiFranco etal., 2011a), but in response 
to hyperpolarizing pulses. It is important to note that 
the kinetic features of di-8-ANNEPS transients associ- 
ated with the activation of IK V reported here are also 
qualitatively different from those reported for the 
other two conductances; whereas the maximal depres- 
sion in those cases were attained almost instantaneously 
after the onset of each hyperpolarizing pulse, the de- 
pression maxima associated with IK V take time to man- 
ifest themselves. As illustrated in the inset to Fig. 7 A, 
after the onset of depolarizing pulses, di-8-ANEPPS 
transients initially raised to a peak, decayed rapidly to 
a minimum as a result of the delayed activation of IK V , 
and subsequendy grew progressively, mainly because 
of the inactivation of IK V . Altogether, the contrasting 
kinetic features between di-8-ANEPPS transients reported 
here, with respect to those associated with the other 
two conductances, are in agreement with differences 
between the time- and voltage-dependent activation of 
gK v with respect to constitutively open conductances. 

Probably the most important finding of the present 
work is precisely that, as illustrated in Fig. 8, the overall 
properties of di-8-ANEPPS transients, including the 
time- and voltage-dependent depression associated with 
the presence of IK V (Fig. 8 A) , can be accurately pre- 
dicted by model simulations (Fig. 8 B) provided that a 
significant fraction of gK v is included in the TTS mem- 
branes (the contrasting case when no channels are 
placed in the TTS is illustrated in Fig. 8 D) . As shown in 
the Appendix, the radial cable model equations, stipu- 
lating the existence of an access resistance (R s ) in series 
with the TTS, were modified from the original equa- 
tions of Adrian and co-workers (Adrian et al., 1969; 
Adrian and Peachey, 1973) to include time- and voltage- 
dependent K v channels (also described in the Appen- 
dix) , similar to what was done before for nonregenerative 
ionic conductances (Heiny et al., 1983; Ashcroft et al., 
1985; DiFranco et al., 2011a) and recently with the Na 
conductance (DiFranco and Vergara, 2011). In the Ap- 
pendix, we explain how the radial cable elements of the 
TTS include K v channels with properties extensively 
described in Fig. 5. In addition, the model contem- 
plates a diffusion equation, modified from Barry and 
Adrian (1973; Friedrich et al., 2001), to calculate 
K concentration changes in the lumen of the T tubules 
resulting from fluxes through K v channels at each 
radial segment of the TTS. Both the radial cable 
model and diffusion equations were integrated si- 
multaneously using conventional numerical methods 
(see Appendix) . 
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Quantitative evaluation of the relationship 
between attenuation and currents 

Similar to our analysis of the attenuation of di-8- 
ANEPPS in the presence of large C1C-1 chloride cur- 
rents (DiFranco et al., 2011a), we developed in this 
paper a method to quantitatively assess the properties 
of the attenuation in di-8-ANEPPS transients recorded 
in the presence of relatively large IK V (Figs. 9 and 10). 
The intrinsic normalization involved in the calculation 
of the percentage attenuation from experimental re- 
cords allowed for quantitative comparisons with model 
predictions of mean TTS potential changes because 
they were normalized following the same definition. 
Thus, our results demonstrating the voltage depen- 
dence of the experimental peak attenuation was not 
only to correlate this property with the activation of 
peak IK V , as illustrated in Fig. 10, but, most impor- 
tantly, to allow for comparison with radial cable model 
predictions of the attenuation, as illustrated in Fig. 11. 
Interestingly, by keeping the IK V predicted by the 
model constant and varying the gK^S/gKrTTS ratios, 
we could examine plausible scenarios that ultimately 
allowed us to quantitatively establish that the distribu- 
tion of gK v in the sarcolemma and TTS membranes 
(gKy-S/gKy-TTS) that is most compatible with the ex- 
perimental data must be within a range between 40:60 
and 60:40. It is remarkable that this range of channel 
densities (per unit surface membrane of each com- 
partment) , including, of course, the likely possibility 
that gK v is equally distributed in the surface and TTS 
membranes, is quite similar to that found previously in 
our laboratory for C1C-1 (DiFranco et al., 2011a) and 
Na v 1.4 (DiFranco and Vergara, 2011). 

Another important inference from the results shown 
in Fig. 11 is that the magnitude of attenuations seen 
experimentally can only result if the contribution 
from the TTS represents a large fraction of the total 
IK V recorded in intact fibers. In fact, for the very 
plausible case analyzed, >70% of the total IK V seem- 
ingly arises from the TTS. It is important to mention 
here that the results from these analyses depend on 
the value of the access resistance to the TTS (R s ) used 
in model simulations. In this study, we used Rs = 40 flcm 2 , 
which is the same used previously regarding the analy- 
sis of C1C-1 distribution (DiFranco et al., 2011a) but is 
significantly smaller than the 120-150 flcm 2 suggested 
for mechanically isolated amphibian fibers (Adrian 
and Peachey, 1973; Kim and Vergara, 1998). How- 
ever, this is slightly larger than the preferred values 
of 10-20 flcm 2 recently reported in our analysis of 
the overshoots associated with the activation of the 
Na conductance (DiFranco and Vergara, 2011). The 
reason for the discrepancy with the latter work is that 
the experiments reported here were all done using 
NMG as the main cation, which affords a reduced 



conductivity (to 0.95 fl/cm) in the external solution, 
likely increasing R s . 

The question that may arise from our results, and 
that we can answer now, is how can the ^56% reduc- 
tion in capacitance afforded by osmotic treatment re- 
sult in ^40% reduction in the mean maximum peak 
IK V ? The answer is given by the model. When gK v is 
equally distributed between the sarcolemma and the 
TTS, the majority of the IK V arises from peripheral 
segments of the TTS. In osmotically treated fibers, 
these TTS peripheral segments that remain connected 
undergo voltage changes more similar to the voltage- 
clamp pulses at the surface membrane and thus con- 
tribute in larger proportions to IK V than those of more 
internal segments of the radial network whose voltage 
changes are significantly attenuated with respect to 
the command pulses. di-8-ANEPPS transients recorded 
from optimally detubulated fibers (unpublished data) 
indicate that voltage changes in peripheral TTS rem- 
nants, together with the sarcolemma, show observable 
deviations from those commanded by the voltage 
clamp. These are significantly less pronounced than 
those observed for the entire TTS network in normal 
fibers but more than what would be expected for the 
sarcolemma itself. This is in contrast with the mem- 
brane capacitance, which, being measured under pas- 
sive conditions, reflects more precisely the proportion 
of the TTS that remains connected. 

K accumulation in the lumen of the TTS 

As explained previously (also see Appendix), our 
model incorporates explicit equations for K diffusion 
in the TTS lumen; this is a highly relevant parameter 
for understanding the kinetics of decay of IK V records 
and the repolarization of the TTS during and after 
long-lasting depolarizations, respectively. We plan to 
take full advantage of this feature of the model in fu- 
ture experimentation and analysis, but for the pur- 
poses of this paper, the changes in the luminal K 
concentration do not play a considerable role. Never- 
theless, as shown in Fig. 8 (A and D) , inward tail cur- 
rents and optical transients slowly returning to 
baseline are clearly observed after the end of the 
pulses. Both of these observations could have re- 
sulted from K accumulation in the lumen of the TTS 
and/or from the presence of unblocked Ca 2+ chan- 
nels. Interestingly, model calculations predict signifi- 
cant [K] accumulation in the lumen of the TTS, 
which, for the largest pulse (to 110 mV) in Fig. 8 D, 
exceeded 27 mM in the inner elements of the TTS; 
however, the simulations paradoxically failed to re- 
produce the slow return to baseline seen experimen- 
tally in optical transients (Fig. 8 A). An obvious 
explanation for this model limitation would be that it 
lacks Ca 2+ conductance. However, the question that 
has yet to be answered is, why doesn't the model report 
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the significant increase in luminal [K]? The reason 
for this stems from the rapid closing of gK v after the 
end of the pulses and the absence of gK IR in the 
model. Namely, without gK JR and gK v , the increase in 
luminal [K] is not being "sensed" by a significant K 
permeability. The aforementioned reasoning consis- 
tently explains the absence of tail currents in model 
simulations (Fig. 8 E) . In contrast, the experimental 
current records (Fig. 8 D) report small but significant 
tail currents that could arise from unblocked Ca 2+ 
channels and/ or incompletely blocked (by 5 mM Rb) 
K R channels; as stated previously, the current model 
does not reproduce either of them. It is important to 
mention that changes in luminal [K] do not signifi- 
cantly affect the analysis followed to determine the 
gK v distribution between the sarcolemma and the 
TTS (as outlined in Fig. 11) because at the peak of IK V 
records, the luminal [K] changed <5 mM for the larg- 
est currents, even in the innermost radial elements of 
the TTS. 

Although several model parameters were adjusted in 
this paper to explain the various features of the experi- 
mental data, the one that needed the most critical at- 
tention to predict simultaneously the properties of IK V 
and the optical data was the relative partitioning of 
gK v between surface and TTS membranes. An impor- 
tant merit of the current work is precisely that we suc- 
ceeded in providing a quantitative assessment of the 
narrow range in the relative distributions of K v chan- 
nels between these two membrane compartments that 
can explain the experimental data. This is a funda- 
mental contribution for the future assessment of the 
electrical properties of skeletal muscle fibers in nor- 
mal and diseased animals. 



APPENDIX 

Model equations for IK V as they are observed 
in osmotically treated fibers 

We assume that the currents recorded from osmoti- 
cally treated fibers arise from surface and peripheral 
TTS membranes that are under voltage-clamp control. 
Thus, IK V records acquired from these preparations 
constitute a reasonable idealization of currents carried 
by skeletal muscle K v channels as expressed in situ. As 
suggested in the body of this paper, we must consider 
the contributions from at least two channel types to 
account for the properties of IK V records. We call them 
channels A and B, and we assume that both of them 
display the general characteristics of Shaker (and/ or 
Sham) K channels; namely, they are multiple-state 
channels that undergo N- and C-type inactivation pro- 
cesses. The analysis of data obtained in osmotically 
treated fiber suggests that channel A undergoes com- 
plete inactivation, whereas channel B behaves as a 



more typical N- and C-inactivating channel. Thus, for 
simplicity, they were modeled as follows. 

Channel A. 

C0 < > C1 < 3 ° A > C2 < 2a " > C3 < "' > Q < *" > IN, 

h <0j Vt 

(SCHEME Al) 

where CO, CI, C2, C3, and O are closed and open states 
of a typical HH K channel (Hodgkin and Huxley, 
1952b) and IN is a terminal inactivated state. The sys- 
tem of differential equations that permit the calculation 
of the contribution of channel A to the total IKy(V,t) 
records are 



dCl 

dt 

dC2 



4a A xC0-p A xCl-3a A xCl + 2P A xC2; 
3a A x CI - 2P A x C2 - 2a A x C2 + 3fi A x C3; 
2a A x C2 - 3P A x C3 - a A x C3 + 4p A x O; 
= a A xC3-4fi A xO-<\> A xO + q> A xIN; (A1) 
= <\> A x O - (p A x IN; 



dl 

dC3 

dt 
dO 

dt 
dl 
dt 

C0 = l-(Cl + C2 + C3 + O + IN); 



S K v\ A =S K v\ A xO;and 
IK v \ A =gK v \ A x(V-E K ). 

The rate constants a and (3 had typical voltage depen- 
dences given by the generic formulas 



a A (V)-. 



a 



' andP A (l/) = P A xe ^ , (A2) 



1- 



V-Va 
Ka 



where €> A and <p A are the open channel (voltage inde- 
pendent) forward and reverse inactivation kinetic rate 
constants, respectively. 

Channel B. Sequential N- and C-type inactivation. 

co < 4a ° > C1 < 3 °° > C2 ^^ C3 < _J!f_> O r JiL ~> /A f <~ J ~^C 

Hi, 2 ft, 3 ft, 4ft,. ~ 

(SCHEME A2) 

The system of differential equations that permit the cal- 
culation of the contribution of channel B to IK v (V,t) 
records are 



DiFranco et al. 131 



= 4a B x CO - P 5 x CI - 3a 5 x CI + 2P B x C2; 

dt 

dCI 

— — — = 3a 5 x CI - 2P fi x C2 - 2a B x C2 + Sp^ xC3; 
rfC3 

= 2a fi x C2 - 3P B x C3 - a B x C3 + 4p B x 0; 

dt 

— = a B x C3 - 4p B x 0 - <|> B x O + (p B x IN; 
dt 

— = <b B xO-q B xIN -yIN + SIC; 
dt 

— = yIN-8IC; 
dt 

C0= l-(Cl + C2 + C3 + 0 + /iV + /C); 



(A3) 



S K v\ B =S K v\ B xO;and 
IK v\b=S K v\b x ( V - E k)' 



where CO, CI, C2, C3, O, and IN are the same as for 
channel A and IC is a C-inactivated state, connected 
with IN through the voltage-independent rate constants 
7 and 8. In this case, the rate constants a and (3 voltage 
dependences are given by 



a B (V)- 



a B (V-V B ) 



and & B (V) = & B x, 



v-v B 



(A4) 



Integration of differential equations 

Following the HH approach for K channel activation, 
the initial conditions (t = 0) were established assuming 
that the open probability at the resting potential (Vh) 
obeys the typical equilibrium equation and that the in- 
activated states are not populated. Thus, 



O(0): 



a. 



a x +p, 



IN(0) = 0 for channels A and B; 7C(0) = 0 for channel B. 
In Eq. A5, ax and (3x are obtained from Eqs. A2 and A4 
for channels A and B, respectively, evaluated at V = Vh- 
From these equations, the probability of all the closed 
states at rest can be calculated as 



C0(0) : 
C2(0) : 



(«x+Px) 4 



, C1(0) : 



4a x P, 



(a x +P ; 



\4 ' 



and C3(0) = 



4a/p, 



Px) 4 



The total IK V records were assumed to be generated by 
the independent contributions of channels A and B; thus, 



IKy=IK v \ A 



IK v \ B ^(gK v \ A+g K v \ B )x(V-E K ) 



(A6) 



Henceforth, the system of simultaneous differential equa- 
tions was numerically integrated using a fourth order 
Runge-Kutta algorithm either in a FORTRAN program 
that calculates the radial cable model predictions or in a 
Berkeley Madonna (Macey and Oster, version 8.3.18, 
2011) program. The parameters a x , V x , k aX , p x , k px , 
<E>x, and 9x (for channels A and B) and 7b and 8 B (for 
channel B) were determined by simultaneously fitting IK V 
records obtained at several voltages from osmotically 
treated fibers with the aid of a Levenberg-Marquardt mini- 
mization algorithm least-square fitting routine included in 
Berkeley Madonna. Mean values are listed in Table Al . 

Radial cable model equations for the TTS of mammalian 
skeletal muscle fibers, including potassium conductance 
and ion diffusion 

The radial cable model equations are essentially the 
same as those described in detail previously (DiFranco 
and Vergara, 2011; DiFranco et al., 2011a), which in 
turn followed the nomenclature published elsewhere 
(Adrian et al., 1969; Ashcroft et al., 1985; Kim and Vergara, 
1998; DiFranco etal., 2007) and the assumption (Adrian 
and Peachey, 1973) that the lumen of the TTS is sepa- 
rated from the extracellular fluid by an access resistance 
(R, in O centimeters squared). The partial differential 
equation that governs the radial- (r) and time ( ^depen- 
dent changes in T tubule membrane potential (u(r,t)) 
in response to voltage changes at the external boundary 
is (Adrian etal., 1969) 

TABLE A1 

Para meter values for Ky channels A and B 



Parameter 



Channel A 



Channel B 



(A5) _ 



(ms" 1 ) 


0.0319 ± 0.009 


0.049 ± 0.0.35 


P* (ms" 1 ) 


0.0142 + 0.004 


0.346 ± 0.043 


V * (mV) 


2.3 ± 1.8 


5.7 ± 3.2 


Kx (mV) 


13.7 ±0.9 


24.2 + 5 


k px (mV) 


14.1 ± 0.9 


33.9 ± 5 


<D X (ms -1 ) 


0.007 ± 0.003 


0.0105 ± 0.003 


cp x (ms -1 ) 


1 x 10 -5 


0.0103 + 0.0016 


7b (ms -1 ) 




0.0021 ± 0.0021 


8 B (ms -1 ) 




1 x 10- 6 ± 1.4 x 10 -5 


S K vl 9 
lx (mS/cm 2 ) 


0.86 ±0.4 


1.92 ±0.4 


Percent contribution 


0.31 ± 0.04 


0.69 ± 0.06 


gieak (mS/cm 2 ) 


0.1 


0.1 



The values correspond to mean ± SEM calculated from fits to records 
obtained for 80-, 100-, 120-, 140-, 160-, and 180-mV pulses from nine fibers. 
Note that the IC state does not exist for channel A. The rate constants a 
and p for channels A and B had voltage dependences given by Eqs. 2 and 
4, respectively. See text for specific deviations from these values. 
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d 2 u 1 du 2 du 

— 5- 4- = v i u + — , 

dR 2 R dR dt 



(A7) 



where a is the radius of the m uscle fiber, R = r/a, T = G L 
t/ ( C w a 2 ) , and v = a JG W / G L . In these equations, the pa- 
rameters C w and G m are the capacitance (in microfarads/ 
centimeters cubed) and conductance (in Siemens/ 
centimeters cubed) of the tubular membrane per unit 
volume of muscle fiber, respectively ( C w = C w p / q and 
G w = G w p / q ) . Also, G L is the effective radial conductiv- 
ity (G L =G L pa, in Siemens/centimeters). Eq. Al must 
be numerically integrated to predict the voltage of TTS 
cable elements when nonlinear conductances, in parallel 
with capacitive elements (Adrian and Peachey, 1973; 
Heiny et al., 1983; Ashcroft et al., 1985; DiFranco and 
Vergara, 201 1 ; DiFranco et al., 201 la) , are assumed to be 
present. To this end, we replace the term xru in Eq. Al 
with the more general term 

2 I w (u,r,t) 

CI , 

G L 

where Iw(u,r,t) is a generalized current normalized per 
volume of fiber. Keeping the same definitions for R and 
T, Eq. Al is transformed into 



d 2 i 



dR 1 



1 du 
R dR 



,2 Iw 



du 
dT' 



(A8) 



For the simulations in this paper, we assumed that 
Iw(u,RT) is contributed by a residual (almost negligi- 
ble) leak current (DiFranco et al., 2011a) but primarily 
by the flow of K ions across the TTS walls through volt- 
age-dependent K channels generally characterized by 
gK v . Thus, 



I w (u, R, T) = H x [I lmk (u) + IK v (u, R,T)], 



(A9) 



where Ii,, a k(u) and Ik{u,R,T) are the leak and K v cur- 
rent per centimeters squared of TTS membrane, re- 
spectively. In this paper, we make the assumption that 
hak(u) = gieak(u - V H ) , where g Uak is constant (see 
Table 1 ) . Also, IK V in Eq. A9 corresponds to the sum 
of the current contributions of channels A and B as 
recorded from osmotically treated fibers and sum- 
marily calculated from Eq. A6. In addition, unless 
otherwise noted, we make the further assumption that 
the proportion of channel A and B contributions is 
the same in the TTS, as it was obtained in detubu- 
lated fibers, an assumption which may or may not be 
fully justified. 



fiber. As described previously (Kim and Vergara, 1998; 
DiFranco et al., 2007, 2011a,b), at a given time j, the 
finite differences approximation of the partial differen- 
tial equation of the T tubular voltage (Eq. A8) , while 
using an implicit Crank-Nicolson algorithm (Crank, 
1975; Gerald, 1978), yields the following equation for 
an arbitrary annulus i: 



v (2i + 1) J /v^n i*xv (2i-1) J +1 
X x x u , - (X + 1) x u +Xx x u,_i = 

4i 4i 



(2i + l) i 
-X x — - x u ] M + (X - 1) x U J . - 



4i 



Xx 



(2i-l) j 2 (IK v )t 



(A10), 



4i 



8T, 



where (IK V )J is the current per unit fiber volume (cal- 
culated from Eq. A9) flowing through the T tubular ele- 
ment at shell i and at the time interval j during the 
numerical integration process. Eq. A10 is a recursive 
formula allowing the calculation of m/ +1 at a time interval 
8 T while knowing uj . The system of tridiagonal coeffi- 
cient matrices was solved using an LU (lower and upper) 
decomposition algorithm (Gerald, 1978). The integra- 
tion of Eq. A6, updating the value of IK v for every cable 
element and for the surface membrane, was performed 
with a fourth order Runge-Kutta algorithm at every 
time step. The values of the general cable parameters 
used for the simulations shown in the paper are sum- 
marized in Table A2. Specific parameters are given in 
the figure legends. 



Potassium diffusion in the TTS lumen 

Changes in luminal potassium concentration in the 
TTS ([K] Q ) occurring in response to the current flow 
across the T tubule walls were calculated from simulta- 
neous integration of the diffusion equation (Barry and 
Adrian, 1973; Friedrich et al., 2001) 



dlK] 0 (r,l)_oD K d 



di 



rd[K] 0 (r,t) 



PF ' 



(All) 



where At is the diffusion coefficient of potassium ions 
in the lumen of the TTS (in centimeters squared/ sec- 
onds) , F is the Faraday constant, Ik is the total potassium 
current (both components added) per centimeter 
cubed of fiber, and p is the fractional volume of the 
TTS. Using the dimensionless variables R = r/a and t = 
<rD K t/a 2 , this equation becomes 



Numerical integration of model equations 

TTS voltage. The TTS cable is assumed to be made of 

n = 60 radial shells, sealed at the center of the muscle 



dR 2 



1 d[K] ri d[K] n 



2 r 

a h 



R dR dx poD A F 



(A12) 
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TABLE A2 



Specific parameter values for passive 

r J r J r 


radial cable 






Syml )< )1 


Value 


Radius (urn) 


a 


20-30 


Specific capacitance TTS wall {(xF/cm 2 ) 


( , vv 


0.9 


XXS lumen conductivity (NMG external) 






(mS/cm) 






TTS access resistance (Ocm~) 


R, 


30-40 


Fraction of fiber volume occupied by the TTS 


P 


0.004 


Volume to surface ratio of the TTS (cm) 


s 


1 x 10~ 6 


Tortuosity factor on the TTS 




0.32 


K diffusion coefficient (cm 2 /s) 


D K 


1 x 10~ 5 



Potassium concentration in the lumen of the TTS. The 
finite difference approximation and Crank-Nicolson 
algorithm used for the integration of Eq. A12 to cal- 
culate the luminal K concentration ([K 0 ]J) at every 
shell i and at time interval j follow an analogous pro- 
cedure to that described for the TTS voltage calcula- 
tions. For simplicity, let's denote [K] 0 = c. Thus, at a 
given time j, the finite differences approximation 
of the partial differential equation for [K] Q in the 
TTS (Eq. A12) is 



BR 2 



1 dc 



— > 



2i8R 



R dR 
{(2i + l) C / +1 



where c\ is the finite difference representation of the 
luminal K concentration at radial position i at time j 
and 8R is the normalized distance between shells. Also, 



8c 
8t 



8t 



-; then 



8t 



8R 2 



,(2i + l) 



-2c/ 



a 2 /, 



(A13) 



poD A -F 



-(2i + l)X x c'H + 4i(X + 1) - (2i - l)X x c>^ 
(2i + 1)X x c{ +1 + 4*(l - X)4 + (2i - 1)X x 



■ 4i x ■ 



a 2 8xI K | 
paD K F 



(A14) 



At the outside boundary (i = n), we apply Fick's law to 
the T tubules' opening: 



/'+! / , 

c = c + 



_5x_ 
5i? 2 



(2n + l)[q], g _,. | (2n-l)cj 1 



2n 



2n 



a 2 8x x I K 
paD K ¥ 



From Crank (1975), for the inside (closed) boundary (i= 0), 



4 



- r ] 

- c o 



45t 
6R 2 



a 8x x I k 
pcD K ¥ 



Eq. A14 is a recursive formula allowing, together with 
the respective equations of the boundary conditions, 
for the calculation of c/ + at a time interval 8 T while 
knowing c{ . The system of tridiagonal coefficient matri- 
ces were solved using an LU decomposition algorithm 
as described for the TTS voltage. 

The values of the general cable and K diffusion pa- 
rameters used for the simulations shown in the paper 
are summarized in Table A2. Specific conductance pa- 
rameters are given in the figure legends. 



Total current and optical signals 

The total ionic currents calculated from the integra- 
tion of the radial cable equations at each radial cable 
element of the TTS was the sum of two contributions: 
(1) a surface membrane component, calculated from 
the direct numerical integration of Eq. A6 in response 
to Vcom; and (2) a TTS component, which encom- 
passes the effective sum of currents originated in this 
membrane compartment; this was calculated from the 
application of Kirchhoff laws at the external opening 
of the TTS (Ashcroft et al., 1985; Kim and Vergara, 
1998) as 



where I K \ J . is the current per unit fiber volume (cal- 
culated through the T tubular element at shell i and 
at the time interval j during the numerical integration 
process). Eqn. A13 allows us to establish recurrence 
using the Crank-Nicolson algorithm. If we denote 



St 
dR 2 



then 



IKy 



R. 



where u{ is the voltage across the outermost segment of 
the TTS at time interval j. 

Under the assumption that di-8-ANEPPS optical sig- 
nals occur homogenously at superficial and deep re- 
gions of the TTS and that changes in its optical properties 
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with the transmembrane voltage occur identically at ev- 
ery submicroscopic region of the TTS, the optical signal 
obtained within our illumination disk is expected to 
represent an ensemble mean of the voltage contributed 
by every cable element in the TTS (AV^) weighted by 
the radius of each annulus. This is calculated, for each 
successive time step, using a numerical trapezoidal inte- 
gration algorithm based on Simpson's rule (Gerald, 
1978) from the formula 




However, as described in detail previously (Kim and 
Vergara, 1998), the actual contributions of different 
radial elements to the overall fluorescence transients 
should be corrected by the optical properties of the ob- 
jective used both to illuminate the preparation and to 
collect fluorescence light. 

Fig. Al shows the normalized fluorescence de- 
tected from a thin film of rhodamine G (~10 jim thick) 
deposited at the bottom of the experimental chamber 
as the z axis distance was varied in 1-um steps to slightly 
above and below the focal plane (z = 0). The continu- 
ous line in Fig. Al was fitted to a Lorentzian function of 
the form 

F-F + 2A ( FWHM } 
" + n Uz 2 +FWHM 2 J 

where F 0 = 0.025, A = 1 19 jim, and the full width at half- 
maximum = 76 um. This function was used to correct 
the weighted mean calculations of the TTS voltage from 
the radial cable model for depth changes in the fluores- 
cence detection. Of the several illumination/detection 
algorithms described by Kim and Vergara (1998), we 




o.o-l , , , 1 , , , , 

-200 -100 0 100 200 

Z distance (pm) 



Figure A1 . Fluorescence intensity depth profile of the objective 
(USAPO 100XO; Olympus) used in the current experiments. 
Data points represent normalized mean fluorescence (three trials). 
The data symbols span the error range. 



used here the global illumination because it is the one 
that accurately represents the experimental conditions 
in the current work. 
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